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ABSTRACT 

Context. The classical T Tauri star (CTTS) AA Tau has presented photometric variability that was attributed to an inner disk warp, 
caused by the interaction between the inner disk and an inclined magnetosphere. Previous studies of the young cluster NGC 2264 
have shown that similar photometric behavior is common among CTTS. 

Aims. The goal of this work is to investigate the main causes of the observed photometric variability of CTTS in NGC 2264 that 
present AA Tau-like light curves, and verify if an inner disk warp could be responsible for their observed variability. 

Methods. In order to understand the mechanism causing these stars’ photometric behavior, we investigate veiling variability in their 
spectra and u - r color variations and estimate parameters of the inner disk warp using an occultation model proposed for AA Tau. 
We also compare infrared Spitzer IRAC and optical CoRoT light curves to analyze the dust responsible for the occultations. 

Results. AA Tau-like variability proved to be transient on a timescale of a few years. We ascribe this variability to stable accretion 
regimes and aperiodic variability to unstable accretion regimes and show that a transition, and even coexistence, between the two is 
common. We find evidence of hot spots associated with occultations, indicating that the occulting structures could be located at the 
base of accretion columns. We find average values of warp maximum height of 0.23 times its radial location, consistent with AA Tau, 
with variations of on average 11% between rotation cycles. We also show that extinction laws in the inner disk indicate the presence 
of grains larger than interstellar grains. 

Conclusions. The inner disk warp scenario is consistent with observations for all but one star with AA Tau-like variability in our 
sample. AA Tau-like systems are fairly common, comprising 14% of CTTS observed in NGC 2264, though this number increases to 
35% among systems of mass 0.7M 0 < M < 2.0M 0 . Assuming random inclinations, we estimate that nearly all systems in this mass 
range likely possess an inner disk warp. We attribute this to a possible change in magnetic field configurations among stars of lower 
mass. 

Key words, accretion, accretion disks - stars: pre-main sequence - techniques: photometric, spectroscopic 


* Based on data from the Spitzer and CoRoT missions, as well as the 
Canada France Hawaii Telescope (CFHT) MegaCam CCD, the Euro¬ 
pean Southern Observatory (ESO) Very Large Telescope, and the U.S. 
Naval Observatory. The CoRoT space mission was developed and op¬ 
erated by the French space agency CNES, with participation of ESA’s 
RSSD and Science Programmes, Austria, Belgium, Brazil, Germany, 


and Spain. MegaCam is a joint project of CFHT and CEA/DAPNIA, 
at the Canada-France-Hawaii Telescope (CFHT), operated by the Na¬ 
tional Research Council (NRC) of Canada, the Institut National des Sci¬ 
ences de l’Univers of the Centre National de la Recherche Scientifique 
(CNRS) of France, and the University of Hawaii. 

** Figures 21-24 are only available in electronic form. 
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1. Introduction 

Magnetospheric accretion models account for most of the ob¬ 
served characteristics of classical T Tauri stars (CTTS), such as 
their strong ultraviolet (UV) and optical excess, as well as broad 
emission lines and forbidden emission lines (Shu et al. 1994; 
Hartmann et al. 1994; Muzerolle et al. 2001; Kurosawa et al. 
2006; Lima et al. 2010). These stars have also been observed 
to have strong magnetic fields (Johns-Krull et al. 2001) and X- 
ray emission (e.g., Preibisch et al. 2005). Magnetohydrodynam¬ 
ics (MHD) models suggest that the stellar magnetosphere trun¬ 
cates the inner circumstellar disk at a distance of a few stellar 
radii, where the disk’s matter stress equals the magnetic stress, 
forcing the material within this truncation radius to be accreted 
onto the star, following the magnetic field lines (Bessolaz et al. 
2008). This forms what are called accretion columns, which are 
responsible for the broad emission lines observed, particularly 
those of the hydrogen Baimer series. When the free-falling ma¬ 
terial reaches the stellar surface, accretion shocks result in hot 
spots near the magnetic poles of the stars, causing the observed 
UV and optical excess. Part of the material from the disk can 
also be ejected in a disk wind, resulting in forbidden emission 
lines. These stars also present a strong infrared (IR) excess with 
regard to the photospheric flux because of the absorption and 
re-emission at longer wavelengths of stellar radiation by the cir¬ 
cumstellar disk. 

Magnetohydrodynamics models (e.g., Kulkarni & Ro¬ 
manova 2009; Romanova et al. 2009; Kurosawa & Romanova 
2013; Romanova et al. 2013) have proposed that CTTS can ac¬ 
crete matter from their circumstellar disks through a stable or an 
unstable regime. In the stable regime, a misalignment between 
the star’s magnetic field and the star’s rotation axis causes the 
material present in the region near the disk’s truncation radius 
to fall onto the star via two main accretion columns, following 
the magnetic field lines connecting the disk to the star (Fig. 1). 
This interaction may also lift dust above the disk midplane, cre¬ 
ating an optically thick warp in the inner disk region near the 
truncation radius, located at the base of these accretion columns. 
At the other end of each accretion column, an accretion shock 
on the stellar surface results in a hot spot, located at a high lat¬ 
itude. Observations of the magnetic fields of CTTS have shown 
that this misalignment between the star’s magnetic field and the 
star’s rotation axis is common (Gregory et al. 2012). 

Magnetohydrodynamics models also predict that the stellar 
magnetic field interacts with the disk at and near the co-rotation 
radius. Differential rotation between the star and the disk, and 
throughout the accretion columns, results in a distortion of the 
magnetic field lines, causing them to inflate and eventually dis¬ 
connect and reconnect, restoring the original field configuration 
(Zanni & Ferreira 2013). 

Accretion through an unstable regime occurs when 
Rayleigh-Taylor instabilities between an accretion disk and a 
stellar magnetosphere cause matter to accrete onto the star via 
streams or tongues in random locations around the star, caus¬ 
ing stochastic photometric variability (see Fig. 1 in Kurosawa & 
Romanova 2013). 

If a star undergoing stable accretion is seen at a medium or 
low inclination, its light curve should show periodic variability 
due to the modulation of its main hot spot. If it is seen at a high 
inclination (yet not so high that the flared outer disk completely 
occults the star), the inner disk warp will occult the stellar pho¬ 
tosphere periodically, causing flux dips in the star’s light curve. 
If a star undergoing unstable accretion is seen at a medium or 
low inclination, its light curve will show stochastic variability 


M a 



Fig. 1: When there is a misalignment between the axis of the 
(predominately dipolar) magnetic field (ji) and the stellar rota¬ 
tion axis (Q), magnetospheric accretion may occur via the stable 
accretion mechanism, where two main accretion columns form 
at the regions where the stellar magnetic pole is closest to the 
disk. Two main hot spots due to accretion shocks are located 
near the magnetic poles. 


due to the many hot spots on the stellar surface (Stauffer et al. 
2014). These hot spots are due to the various accretion shocks 
around the star. If these stars are seen at high inclinations, we 
can expect to observe stochastic occultations of the photosphere 
by dust lifted above the disk plane near the base of the accretion 
tongues whenever one passes in front of our line of sight, leading 
to light curves dominated by aperiodic extinction events. These 
aperiodic extinction events could also be due to occultation by 
dust that is lifted above the disk plane as a result of magnetoro- 
tational instabilities (Turner et al. 2010). 

During a study of the photometric and spectroscopic vari¬ 
ability of the classical T Tauri star AA Tauri, conducted between 
1995 and 2007 (Bouvier et al. 1999, 2003, 2007), a simple ge¬ 
ometrical model was presented to explain the behavior of AA 
Tau’s light curve, characterized by a relatively constant level of 
brightness interrupted by quasi-periodic dips of around 1.4 mag¬ 
nitudes (see Fig. 9 in Bouvier et al. 1999). This variability was 
attributed to periodic occultations by an inner disk warp spatially 
associated with two major accretion columns, as in the stable 
accretion regime. This result provided great support for MHD 
models. 

It is interesting to note that after over 20 years of observa¬ 
tions often showing the same quasi-periodic behavior, A A Tau’s 
light curve suddenly changed quite drastically when its average 
brightness level decreased by over 2 magnitudes in 2011 (Bou¬ 
vier et al. 2013), at which point it lost its periodicity. This in¬ 
dicates that, though the inner disk warp structure may be stable 
over the course of many years, it can disappear in a timescale of 
less than one year. 

There are many possible scenarios in which a warp can ap¬ 
pear within the circumstellar disk of a T Tauri star besides the 
one described here (see, e.g., Terquem & Papaloizou 2000; Fla¬ 
herty & Muzerolle 2010). In this paper, we treat the case of clas¬ 
sical T Tauri stars, which are actively subject to ongoing mag¬ 
netospheric accretion, with inner disk warps located at or near 
the disk’s co-rotation radius (see discussion in Sect. 3.5). It is 
reasonable to assume that a warp at this location is due to the in¬ 
teraction between an inclined magnetosphere and the inner disk, 
and therefore associated with the accretion columns, as was con¬ 
sidered for A A Tau in Bouvier et al. (2007). For the stars in this 
paper, we will consider this to be the main mechanism responsi¬ 
ble for the inner disk warp. 
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1.1. The occultation model proposed for AA Tau 

Bouvier et al. (1999) considered that the dynamic interaction be¬ 
tween the inclined magnetic field of the star and its inner disk 
region would cause the material in this part of the disk to be 
lifted, preferentially where the stellar magnetic dipole is closest 
to the disk, resulting in two arch-shaped walls located opposite 
to each other. This warp, located at or near the co-rotation radius 
(R co ) of the accretion disk, eclipses the stellar photosphere when 
it passes through the observer’s line of sight. It has an azimuthal 
extension of 0 W and a maximum height of h max from the disk’s 
mid plane, which occults the stellar photosphere at phase 0.5, 
provoking a minimum in the light curve. 

The model uses simple geometric principles to reproduce a 
synthetic light curve where an optically thick clump of given 
maximum height h max and azimuthal extension 0 W , present at 
the co-rotation radius R co , occults the stellar photosphere of a 
system seen at inclination i. The height of the clump is assumed 
to vary according to 


^( 0 ) — hmax 


COS 


/ 7T(0 — 0p) \ 

l K ) 


(l) 


for values of -0 w /2 < 0 < 0 w /2, where 0o is the azimuth of 
maximum disk height. The clump height is equal to zero for all 
values of |0| > 0 w /2. In the case of AA Tau, the warp found to 
best reproduce its light curve had a maximum height of h max = 
0.3 R co , and an azimuthal extension of 360°. 


1.2. AA Tau-like light curves in NGC 2264 

In a study of CTTS in the young cluster NGC 2264 using light 
curves from the CoRoT satellite (Convection, Rotation and plan¬ 
etary Transits, Baglin et al. 2006), various CTTS with light 
curves similar to that of A A Tau were found (Alencar et al. 
2010). NGC 2264 is a young stellar cluster of approximately 3 
Myr situated in the Orion spiral arm, around 760 pc from the Sun 
(for a review of the cluster, see Dahm 2008). Of the stars in the 
region observed by the CoRoT satellite in 2008, 301 were con¬ 
firmed as members of the cluster, and of these 83 were classified 
as CTTS by Alencar et al. (2010). 

The 83 CTTS were separated into three groups, according to 
the morphology of their light curves. The first group consisted 
of 28 periodic, nearly sinusoidal light curves, whose variabilities 
were attributed to long-lived cold spots; the second group con¬ 
sisted of 23 light curves with the same characteristics as those of 
A A Tau, periodic but with varying amplitude and shape of dips; 
and the third group consisted of 32 non-periodic light curves. 
Two examples of the AA Tau-like light curves, shown in full and 
folded in phase, can be seen in Fig. 2. 

In the same study that identified these light curves, Spitzer 
IRAC single epoch photometry was used to determine whether 
there is dust in the inner disk region of these stars, by identi¬ 
fying near-infrared excess emission. Of the 83 CTTS, 68 were 
also observed by Spitzer. All of the stars with AA Tau-like light 
curves that were observed by Spitzer and CoRoT were shown to 
present dust in the inner disk region. This is consistent with the 
assumption that the photometric behavior of these stars is due to 
obscuration by optically thick material in their inner disks. 

In order to better understand the physical processes that gov¬ 
ern CTTS variability, a new observational campaign of NGC 
2264 was organized in December 2011. The Coordinated Syn- 
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Fig. 2: Two examples of A A Tau-like light curves found for 
CTTS in NGC 2264 by Alencar et al. (2010). (c) and (d) show 
the light curves (a) and (b), respectively, folded in phase. Differ¬ 
ent colors correspond to different rotation cycles. 


optic Investigation of NGC 2264 (CSI 2264, Cody et al. 2014) 1 
is composed of simultaneous observations on fifteen telescopes 
(eleven ground-based and four space-based telescopes), includ¬ 
ing the CoRoT satellite, the Spitzer Infrared Array Camera 
(IRAC, Fazio et al. 2004), the Canada-France-Hawaii Telescope 
(CFHT) MegaCam (Boulade et al. 2003), and the FLAMES 
multi-object spectrograph on the Very Large Telescope (VLT, 
Pasquini et al. 2002). This campaign provided us with simultane¬ 
ous photometric and spectroscopic information in a wide range 
of wavelengths. New CTTS were identified according to their 
UV excess (Venuti et al. 2014), and various new stars with AA 
Tau-like light curves were found, in addition to stars with newly 
identified types of photometric behavior (Stauffer et al. 2014, 
2015). 

Because of the large number of previous studies of NGC 
2264, it was possible to establish reliable criteria for cluster 
membership, as discussed in Alencar et al. (2010). As probable 
cluster members, we used photometric H a and variability from 
Lamm et al. (2004) following their criteria; X-ray detection from 
Ramirez et al. (2004) and Flaccomio et al. (2006) and location 
on the cluster sequence in the (/, R -1) diagram, when available; 
spectroscopic H a equivalent width greater than a spectral type 
dependent threshold (White & Basri 2003); and H a emission 
line width at 10% intensity greater than 270 km/s, as proposed 
by White & Basri (2003) to identify accreting T Tauri stars. 

In this study we present a multi-wavelength analysis of the 
CTTS members of NGC 2264 that appear to have extinction 
dominated photometric behavior. In Sect. 2 we present the obser¬ 
vations used in this study, while in Sect. 3 we show the analysis 
of our data and our results. In Sect. 4 we present a discussion on 
the proposed cause of these stars’ photometric variabilities, and 
in Sect. 5 we present our conclusions. 


1 The CoRoT and Spitzer light curves for all probable NGC 2264 
members, as well as our broad band photometry for these stars, 
are available at http: //irsa. ipac. caltech. edu/data/SPITZER/ 
CSI2264. 
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2. Observations 

2.1. CoRoT data 

The first CoRoT short run occurred in March 2008, when the 
young cluster NGC 2264 was observed for 23 days uninterrupt¬ 
edly. Nearly four years later, during the CSI 2264 campaign, the 
cluster was observed once again by CoRoT from 2011 Decem¬ 
ber 1, to 2012 January 3, during the satellite’s fifth short run. 
CoRoT’s field of view is larger than the full extent of the cluster, 
which fits entirely into one of the CCDs designed for studying 
exoplanets. Stars of magnitude down to R ~ 18 were observed, 
with a cadence of 512 seconds. In some cases the data were taken 
with exposure times of 32 seconds to provide high cadence light 
curves. All of the light curves from the first epoch were re-binned 
to 512 seconds. 

CoRoT observes in a white light bandpass of approximately 
3700A to 10000A. It can, in principle, provide three-color pho¬ 
tometry by passing the stellar light through a low-resolution 
spectral dispersing prism, but for the purpose of this paper we 
used only the integrated white light flux. 

After undergoing standard CoRoT pipeline procedures, de¬ 
tailed in Samadi et al. (2006), the data were subsequently pro¬ 
cessed to remove unwanted effects, such as those caused by 
Earth eclipses, and to remove measurements affected by the 
South Atlantic Anomaly, using a sigma-clipping filter. Care was 
taken not to remove flaring events. In some cases, jumps due to 
detector temperature jumps were present in the photometry, and 
when possible they were removed manually (for more details on 
the correction of CoRoT light curve systematics, see Cody et al. 
2014). 

2.2. Spitzer IRAC data 

The Spitzer Space Telescope mapped a region of ~ 0.8° x 0.8° 
centered at R.A. 06:40:45.0, declination +09:40:40, from 2011 
December 3, to 2012 January 1. Targets were observed approxi¬ 
mately twelve times a day in the IRAC 3.6pm and 4.5yum chan¬ 
nels in its Warm Mission mode. Most of the objects have data 
from both bands, but around 40% of the objects fall near the 
edges of the mapping region, and only have data from one band. 
A staring mode was also used around the beginning of the run 
for a region near the center of the cluster. A cadence of approx¬ 
imately 15 seconds was obtained during four blocks of 20, 26, 
16, and 19 hours, on 2011 December 3, December 5-6, Decem¬ 
ber 7-8, and December 8-9. 

The IRAC data reduction procedures and production of 
IRAC light curves are explained in detail in Cody et al. (2014) 
and Rebull et al. (2014). 

2.3. CFHT MegaCam data 

NGC 2264 was observed using the Canada France Hawaii Tele¬ 
scope’s wide-held imager MegaCam (Boulade et al. 2003) be¬ 
tween 2012 February 14 and 2012 February 28. The MegaCam’s 
held of view spans a full lxl square degree, covering the cluster 
with only one pointing. Data were taken in the u and r bands on 
30 epochs. Reduction of these data is explained in Venuti et al. 
(2014). 

2.4. USNOdata 

On the order of 900 epochs of I-band photometric observations 
of NGC 2264 were obtained at the US Naval Observatory 40" 


telescope from December 2011 to March 2012. A held of view 
of 23’x23’, centered approximately on the cluster center, was ob¬ 
served. The data underwent standard reduction procedures using 
bias and dome hat images that were obtained at the beginning 
of each night. Aperture photometry was performed for all stars 
in each image, and a set of non-variable stars was identihed in 
order to establish zero-points for each CCD frame and construct 
light curves of the targets using differential photometry. 

2.5. Spectroscopy 

For some stars in our sample, 20-22 epochs of VFT FFAMES 
spectra were obtained between 2011 December 4, and 2012 
February 29. Of these, we have four to six epochs of simulta¬ 
neous FFAMES spectroscopy and CoRoT photometry. 

FFAMES can obtain medium- and high-resolution spectra 
for multiple objects over a held of view 25 arcmin in diame¬ 
ter. Since this is much smaller than the extent of NGC 2264, 
even with two pointings, only the central part of the cluster 
was observed (one pointing centered at R.A. 06:41:04.8, dec¬ 
lination +09:45:00, the other centered at R.A. 06:40:58.8, dec¬ 
lination +09:21:54). We used the HR15N grating of the GI¬ 
RAFFE/MEDUSA mode, which covers H a and Fi 6707A with 
a resolution of 17,000. A few of the brighter stars were observed 
with the UVES (Ultraviolet and Visual Echelle Spectrograph) 
red mode, which covers the region from 4800A to 6800A with a 
resolution of 38,700. Only one of the stars treated in this paper 
was observed with UVES. 

Where FFAMES spectra were not available, we used high- 
resolution Hectochelle spectra from 2004-2005 provided by Ga¬ 
bor Furesz (Furesz et al. 2006). They used the 190A wide order 
centered at H a, with a resolution of R~ 34000. These spectra 
were used for eight stars. 

3. Data analysis and results 

3.1. Initial selection of candidates for extinction dominated 
light curves 

Throughout this paper we are concerned with stars that present 
extinction events due to circumstellar material eclipsing the stel¬ 
lar photosphere. This requires that the stars in our sample be 
classical T Tauri stars, with at least some optically thick dust 
still present in the disk. Therefore, after selecting the probable 
members of NGC 2264 following criteria described in Alencar 
et al. (2010), a subset of stars with CoRoT observations in either 
epoch was classified as CTTS based on the following criteria: 
H a width at 10% intensity greater than 270 km/s, H a equivalent 
width higher than a spectral type dependent criterion proposed 
by White & Basri (2003), or U-V excess greater than a spectral 
type dependent threshold described in Alencar et al. (2010). A 
new group of CTTS candidates proposed by Venuti et al. (2014) 
based on UV excess measured with the CFHT Megacam during 
the CSI 2264 campaign was also taken into consideration. We 
also consider the two stars Mon-56 and Mon-14132 as candi¬ 
date CTTS since they have IR excess and photometric variabil¬ 
ity characteristic of the CTTS studied in this paper, though we do 
not have enough spectral information to classify them as CTTS 
and they present no considerable UV excess. 

Of the ~ 300 confirmed members observed in 2008, we clas¬ 
sified 95 as CTTS or CTTS candidates, 12 more than had been 
considered by Alencar et al. (2010). Of the - 500 confirmed 
members observed in 2011, 148 were classified as CTTS or can¬ 
didates. Of all of these CTTS and CTTS candidates, 84 were ob- 
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served in both epochs. Therefore our initial sample of possible 
CTTS is of 159 cluster members, observed with CoRoT in either 
2008, 2011, or both. Sousa et al. (in prep.) performed an in-depth 
study of all of these CTTS, focusing on the VLT spectroscopy. 

The CoRoT light curves of these 159 stars were separated 
into different groups based on their visual morphologies. This 
was done based solely on the CoRoT light curves, with no color 
information. One of these groups consists of AA Tau-like light 
curves, which present a relatively stable maximum brightness 
level interrupted by quasi-periodic flux dips that are neither spot¬ 
like nor eclipse-like in nature. The term quasi-periodic means 
that the flux dips present a stable period but vary in shape and/or 
amplitude between rotation cycles. They vary by on average 25% 
in width and 30% in amplitude from one rotation cycle to the 
next. Their periods range from 3 to 10 days and are consistent 
with Keplerian rotation near the co-rotation radius (see Fig. 13 
and discussion in Sect. 3.5). They are usually accompanied by 
other narrower, shallower flux dips that are not necessarily pe¬ 
riodic and can have significant structures within the main flux 
dips. 

Bouvier et al. (2007) showed that the photometric behavior 
of the star AA Tau varied significantly between observations sep¬ 
arated by a few years, showing at times the quasi-periodic flux 
dips attributed to the inner disk warp, at times more than one 
flux dip per stellar rotation period, and at times with almost no 
variability (see Fig. 15 in Bouvier et al. 2007). The light curves 
we have described as AA Tau-like are based on the observations 
of AA Tau from 1995. There are, however, some light curves 
that are similar to the 1995 AA Tau light curve, but that either 
do not have a well-defined maximum brightness, or show rela¬ 
tively weak periodicity (for example, see star Mon-1054 in Fig. 
20). Cody et al. (2014) do not classify these as AA Tau-like light 
curves, but they are important in the context of this paper because 
they may represent systems that are undergoing events similar to 
those that AA Tau underwent during its other phases. Therefore, 
we consider stars that present these light curves in one CoRoT 
observation or the other to be candidate AA Tau systems along 
with the stars that have classical AA Tau-like light curves. We 
analyze their photometric variability in the same way as we do 
the A A Tau-like light curves, in order to establish whether the 
stable accretion mechanism viewed at a high inclination could 
also be responsible for their variability. 

Another group that is of interest to this study is composed of 
light curves with similar characteristics to the AA Tau-like light 
curves, but that show no obvious periodicity. These light curves 
appear to be dominated by aperiodic extinction events. These 
two groups can easily be distinguished from the periodic, nearly 
sinusoidal, variability attributed to long-lived cold spots, or the 
nearly symmetric, narrow flux bursts described in Stauffer et al. 
(2014) and attributed to stochastic accretion events. A few ex¬ 
amples of different types of variability are shown in Fig. 3. For 
a detailed discussion of these and other light curve morpholo¬ 
gies, including a statistical analysis of symmetry in the flux and 
periodicity, see Cody et al. (2014). 

Fig. 4 shows the CoRoT light curve and USNO I-band light 
curve, shifted to coincide in magnitude, of two AA Tau-like 
stars. It is evident from this figure, though the I-band photome¬ 
try coincides closely with the CoRoT photometry, that the lower 
cadence makes it difficult for AA Tau-like behavior to be iden¬ 
tified using only ground-based data. This shows the importance 
of high cadence space-based photometry to better understand the 
phenomena that occur among these young stellar objects. 

Table 1 shows information on the stars in our sample that are 
considered to have AA Tau-like or aperiodic extinction domi¬ 


nated light curves in either epoch. The CSIMon ID is an internal 
naming scheme devised for the CSI 2264 campaign (see Cody 
et al. 2014), and comprises all cluster members, candidates and 
field stars in the NGC 2264 field of view (for brevity, throughout 
the text, we eliminate ‘CSI’ and the leading zeros and replace, 
for instance, CSIMon-000250 with Mon-250). The CoRoT ID is 
the identification number assigned by the CoRoT satellite to all 
objects it has observed. In some cases where the same star was 
observed in both epochs, two different CoRoT IDs were assigned 
to the same star. For these stars we present the ID from the 2008 
observation. The 2MASS identification is also presented. 

The sixth and seventh columns of Table 1 show how the light 
curve was classified in one epoch or the other. Those labeled 
with AAT present AA Tau-like light curves, those labeled with 
AATc are those that were not classified as AA Tau-like by Cody 
et al. (2014), but that we consider to be candidate A A Tau sys¬ 
tems, and those labeled as Ap show aperiodic light curves. For 
six stars, a possible AA Tau-like variability was found in both 
epochs, meaning the mechanism responsible for this variability 
can be stable on a timescale of a few years. Four stars that present 
AA Tau-like variability were only observed in one epoch. There 
are thirteen stars that were observed in both epochs and classified 
as AA Tau-like in one and aperiodic in the other. Of these, five 
displayed A A Tau-like behavior in 2008, yet in 2011 presented 
aperiodic variability, while the other eight showed the opposite 
trend. This points to a possible change from a stable accretion 
regime to one of unstable accretion (or vice-versa) in a matter 
of less than four years. This is discussed in Sect. 3.3. Another 
five stars were classified as having an aperiodic extinction dom¬ 
inated light curve in both epochs. Five more stars with aperiodic 
extinction dominated light curves were only observed in 2011. 

It is interesting to note that, though many stars suffered a 
transition between AA Tau-like and aperiodic extinction dom¬ 
inated light curves, no star whose light curve was classified as 
AA Tau-like in one CoRoT observing run was classified as spot¬ 
like in the other. This could be due to inclinations, since high 
inclinations would favor extinction events and possibly hide any 
variability due to spots on the stellar surface. 

Another group of stars with what appear to be short duration 
extinction events present in their CoRoT light curves is studied 
by Stauffer et al. (2015). These stars present narrow flux dips 
that are generally shallower than in A A Tau-like light curves, 
and nearly Gaussian shaped, while AA Tau-like light curves 
usually show very irregularly shaped dips, often with signifi¬ 
cant structure in the minima. The narrow dip stars are periodic 
over timescales of tens of days, with periods that are consistent 
with Keplerian rotation periods of the inner disk region. Stauf¬ 
fer et al. (2015) argue that the inner disk warp is not adequate 
to explain these flux dips because of their short durations and 
regular shapes, and propose alternate scenarios to account for 
the variabilities. Two of the stars in our sample of A A Tau-like 
light curves are considered in their paper as well. They are Mon- 
56 and Mon-1131, and they present a rather complex photomet¬ 
ric behavior. Therefore they are studied therein through different 
points of view, in order to investigate other possible causes for 
their variabilities. There are also some stars with A A Tau-like 
light curves that at times present these narrow flux dips along 
with the broader ones. These features of their light curves are 
also discussed in Stauffer et al. (2015). 


Article number, page 5 of 27 


A&A proofs: manuscript no. AATau-like_in_NGC2264 


ft) 

> 

5-h 

G 

o 

x 

w 

73 

O 

G 

a 

S 

o 

73 

G 

.2 

o 

.2 

s 


■8 


G 

> 


2 

C/3 

S' 

C/3 

3 

PlT 


.2 3 

05 G 

?> & 


05 

73 o 

W Cxi 

Ph w 


G O 

73 o 

W Cx| 

Ph W 


ft) _ 
GhP^ 

u ° 

rj w 


CD _ 

§o 

u 8 

P w 


ft) 

Cft 

>n 

Oh 

CO 


04 04 
+1 +1 
O 


<<<<<<<<*><<< 


i vd vd r- 
i P P O 


CO 04 
P 04 


qina;^ 
P 04 04 P 


q in q 

CO i-H 1 -H 


h q 

04 * CO 


, CO 


O O Tl" OO [Q 04 
o cs h o s ' " 


o~ 

d 


i i 


d-Od-^l-^HOoin^Hm-H(N-H 

q q q oj vd rn if) q q q vq q 

1 -H 04 04 04 04 1 -H i-H i-H (N i-H 1 -H i-H 


04 04 Tl" CO »n OO 04 

OV O' in Tl" 04 O O 

O P O P t-h i-h 04 


ON CO O oo O' CO O' O O' 04 O 

(NhncoroM)0'sj-ooo4in 


vd oo q 

04 i-H 04 
+1 +1 +1 
CO OO O 

P d P 

H CO 04 


rv CO H 

o ^ ^ I 
-; oo co 


04 p. 




+i 

o 


VD oo 
OO VO 


in q 

i-H 04 
+1 +1 
cq i-H 

av P 

^h 04 


q in q 

CO 04 04 
+1 +1 +1 
q q q 

00 H 

H 04 04 


cq q q oo q 

i-H i-H i-H 04 CO 

+1 +1 +1 +1 +1 

q in oo i-H q 

d cb P d oi 

CO rH rH 04 CO 


’ P 

04 


co o i-H vp , 

04 CO 04 
+1 +1 +1 +1 , | 
q o q o J 

ON 7t" 1-H i-H 

^H 04 ^H -H ^ 


I I 


I I 


I I I 


00 

P 


I I 


OO OO O 
O H QV 


oo in oo 
I o o q | 
oo On VD 


I I I 


VD 

CO 

vd in 

VD 

in 04 

q 

1 p 

1 q 04 | 

1 3 

| O OV 

P 


P in 

CO 

P OV 


.h H 




< < 


p ft p 

c < < * 


I I I 


eft Gh GhH ^ d 

<<<< 

< 


%. m VD 

£ Itl I 

P VD VD 


. ah H H 

•^5 S3 1 


r U ftH 

H ft H ft 

i< 3 < 


Gh Gh 

< < 


% 

< 


Gh H 

< ^ 


Gh 

< 


3 -^ 


Gh 

< 


M ^ l—M M ^ L J 








_ ^ _ OOI04TtiT)OhO> 


Hh 

c 3 

d 


q 

t _ ( 


CO 

7 t 

04 

04 

VD 


04 

cn 

04 

O) 

q 

q 

OV 

OO 

CO 

CO 

04 


q 

m 

CO 

VD 

CN 

O' 

CO 

q 

q 

in 

q 

O 

Cn 

Eo 

ft) 

G 

ft) D? 

G 1 -H 
.2 2 

o 

04 

> c 3 ^ -G 

2 ft> x 

P 5 _i O ft) 



P 

d 

CO 

p 

P 

vb 

d 

P 

vb 

P 

cn 

P 

oo 

cn 

q 

in 

oo 

p 

p 

P 

q 

CO 

d 

vb 

in 

1 —H 

Ov 

CN 

P 

CO 

P 

oi 

vb 

u 

q 

05 

G 

73 

q h 

'G 

PP 

M O G X 
c C^ 0 X 

G 73 G P 

G d § SP 

q -2 d g 

G 

U 

w 

7 )- 

d 

04 

q 

in 

in 

P 

7 t 

oo 

04 

04 

04 

in 

04 

P 

7 j- 

d 

>n 

P 

in 

vb 

04 

cn 

04 

bv 

P 

in 

bv 

bv 

CO 

oo 

7 j- 

vb 

m 

o 

o 

P 

q 

oo 

q 

p 

q 

P 

q 

d 

7 )- 

o 

in 

P 

vb 

q 

P 

in 

P 

04 

o 

CN 

P 

in 

/-H 

vb 

q 

r- 

2 « 
T Oh 

G 

05 

cP 

Q 

>n 

7 f 

7 t 

in 

04 

cn 

in 

Tt; 

in 

CO 

CO 

7 t 

q 

co 

04 

04 

CO 

CN 

co 

q 

CO 

co 

co 

q 

>n 

in 

q 

04 

co 

q 

04 


hH 

C 

ft) 

G 

G “ 

G 

H 


6 v 

6 v 

6 v 

6 v 

6 v 

6 v 

6 v 

6 v 

6 v 

6 v 

6 v 

6 v 

6 v 

bv 

bv 

bv 

bv 

bv 

bv 

bv 

d 

d 

d 

d 


d 

bv 

bv 

bv 

bv 

bv 

o 

bv 

G 

w) .b 

O 

< 

< 


+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

_G 

5 h 

G 

H 

o 

W) ft) 
G G 3 
73 •h* 

rH M 

12 o 
2 

G 05 

G 

£ 

I) 

° 3 73 q 
x p o q . 
c o G oo 

-g 


VD 

7 j- 

o 

OV 

VD 

OV 

7 t 

CO 

O' 


o 

7 t 

o 

m 

o 


oo 


7 j- 


7 t 

04 

VD 

04 

04 


O' 

OV 

O' 

in 

7 t 

VD 

m 

O 

o 

U 

05 

X) 

ft) G Uh O' 

S gins 

73 •§ o u ^ 


C' 

oi 

co 

oo 

p 

m 

o 

in 

d 

in 

in 

P 

CO 

OV 

O' 

vb 

q 

P 

q 

oo 

in 

in 

P 

iq 

P 

in 

OV 

in 

o 

q 

oo 

in 

00 

oo 

CN 

cn 

00 

vq 

in 

OO 

o 

vq 

CN 

m 

vb 

q 

P 

q 

o 

cn 

oo 

p 

q 

cn 

04 

Ov 

OV 

OV 

>n 

cn 

bJD 

o 

rn 

C /3 

•n 

in 

o 

q 

in 

o 

7 j- 

o 

7 f 

o 

m 

o 

q 

o 

q 

in 

o 

m 

o 

7 f 

o 


in 

o 

q 

CN 

q. 

m 

o 

7 f 

o 

q 

CN 

CO 

o 

04 

o 

CO 

o 

CO 

bv 

q 

bv 

04 

o 

CO 

o 

q 

bv 

>n 

bv 

q 

bv 

q 

bv 

q 

G 

G 

ft) 

u 

G 2 

d 2 


£ ' o S l § 

c 3 

HH 

tj- 

7 j- 

7 j- 

7 j- 


7 j- 


7 j- 

7 t 

7 j- 

7 j- 

t)" 

7 j- 

7 t 

7 j- 

7 t 

7 j- 

7 j- 

t)" 

7 j- 

7 t 

7 j- 

7 f 

7 j- 

7 j- 

CO 

CO 

7 j- 

7 t- 

CO 

CO 

CO 

CO 

2 

O G 

S 

OO 

o 

o 

04 

H S. 3 - ” o 

cn 

CD 


vb 

7 )- 

vb 

vb 

vb 

CO 

vb 

7 j- 

vb 

vb 

vb 

in 

vb 

o 

vb 

04 

vb 

CO 

vb 

04 

vb 

OV 

vb 

7 j- 

vb 

O' 

vb 

OV 

vb 

oo 

vb 

CO 

vb 

CO 

vb 

04 

vb 

vb 

OV 

vb 

m 

vb 

04 

vb 

VD 

vb 

04 

vb 

VD 

vb 

CO 

vb 

o 

vb 

7 |- 

vb 

m 

vb 

O' 

vb 

ov 

C /3 

-G 

cS 

73 

ft) 

> 

& 8 

. 25 m a £5 

x ^ q q ^ 

ft) ft) ^ ft) » 

O X ^ QQ O 




o 

CO 

O' 

7 t- 

VD 

o 

O' 

VD 

04 

CO 

7 j- 

oo 

CO 

i-H 

m 

oo 

7 j- 

O' 

O' 

i-H 

04 

o 

VD 

m 

i-H 

OV 

04 

7 |- 

CO 

7 t 

04 

m 

ft) 

$-H 

PLh X 


2 



04 

o 

in 

7 j- 

04 

04 

04 

tT 

m 

in 

04 

04 


in 


CO 

7 |- 

in 

o 

7 t 

7 j- 

tJ- 

7 |- 

t|- 

wo 


7 j- 

m 

04 

04 



<+H 

ft) 

q s 


r G 


o 

i-H 

m 

7 j- 

OO 

04 

m 

O' 

o 

7 j- 

VD 

CO 

OV 

7 t- 

OV 

OV 

oo 

VD 

o 

O' 

OO 

7 j- 

7 |- 

OV 

o 

7 f 

VD 

0 ' 

7 t 

O 

O' 

in 

VD 

o 

C /3 


G 

H 


m 

7 j- 

7 t 

m 

04 

CO 

m 

7 j- 

m 

co 

CO 

7 j- 

7 j- 

co 

04 

04 

CO 

04 

CO 

7 |- 

CO 

CO 

CO 

7 |- 

m 

m 

tJ- 

04 

co 

7 |- 

04 

o 

i-H 


G 

'S (L) 

g q 

Q 

q 2 *h ^ 


cn 

OV 

OV 

OV 

OV 

OV 

OV 

OV 

On 

OV 

OV 

OV 

OV 

OV 

OV 

OV 

OV 

OV 

OV 

OV 

OV 

OV 

OV 

OV 

On 

OV 

OV 

OV 

OV 

OV 

Ov 

OV 

o 

OV 

G 

ft) 

Q 

H—\ -2 ^ ‘H H-H 


cn 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

O 

o 


o 

O 

05 

O 73 

H—H 

ft) 

O 05 1-H 05 , 


< 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

G 

G 

C 2 g s % g 


s 

m 

CO 

o 

OV 

VD 

OV 

7 j- 

CO 

O' 

i—H 

OV 

7 t 

o 

m 

OV 

i—H 

oo 

i—H 

7 j- 

1 —H 

CO 

04 

VD 

04 

04 

1 —H 

O' 

OV 

O' 

7 j- 

7 j- 

VD 

m 

o 


<Z? § 

rG 

ft) 

G 

-h 6 M u n 


i—H 

OO 

m 

in 

i—H 

in 

co 

O' 

OV 

i—H 

o 

m 

in 

O' 

7 j- 

m 

O' 

i—H 

oo 

04 

oo 

VD 

OO 

VD 

m 

7 j- 

7 j- 

o 

oo 

O' 

04 

Ov 

m 

Gh 

Vh 

Hs-H 05 

73 pH ft) ft) 

2 P 2 o & 


04 

CO 

O' 

o 

o 

m 

7 f 

OV 

VD 

7 t 

i-H 

oo 


O' 

7 j- 

OV 

o 

oo 

m 

oo 

CO 


in 

o 

04 

VD 

7 |- 


CO 

O' 

CO 

OV 

Ov 

co 


G 

ft) 



m 

in 

o 

in 

7 |- 

7 f 

m 

o 

o 

o 

in 

m 

7 j- 


in 

o 

04 

i—H 

m 

7 j- 

o 

04 

CO 

04 

CO 

co 

7 |- 

04 

CO 

o 

m 

O 

o 


05 

O 

8 § g c ^ 3 



i—H 

o 


o 

o 

o 

o 




o 

o 

o 

i—H 

o 


i—H 

i—H 

o 

o 


i—H 

o 

O 

o 

OV 

OV 

o 

o 

OV 

OV 

OV 

OV 

ft) 

G 

p .2 

_G 



7 t 

7 |- 

7 t- 

7 |- 

7 j- 

7 f 

7 j- 

7 j- 

7 )- 

7 )- 

7 |- 

7 j- 

7 |- 

7 |- 

7 |- 

7 l- 

7 t 

7 |- 

7 |- 

7 |- 

7 l- 

7 j- 

7 |- 

7 |- 

7 t- 

CO 

CO 

7 |- 

tJ- 

CO 

CO 

CO 

CO 

73 

G 




VD 

VD 

VD 

VD 

VD 

VD 

VD 

VD 

VD 

VD 

VD 

VD 

VD 

VD 

VD 

VD 

VD 

VD 

VD 

VD 

VD 

VD 

VD 

VD 

VD 

VD 

VD 

VD 

VD 

VD 

VD 

VD 

VD 


ft) 




Q 

O 

O 

O 

VD 

O 

OO 

O 

OV 

O 

O' 

O 

O 

O 

O 

OV 

O 

O 

O 

O 

OO 

O 

m 

O 

O 

O 

OV 

O 

O 

O 

CO 

o 

o- 

o 

O' 

O 

7 j- 

O 

O 

o 

CO 

O 

OO 

O 

O 

O 

OV 

O 

OO 

O 

in 

o 

oo 

O 

O 

OO 

O 

O' 

o 

O' 

O 

O 

o 

7 l- 

ft) 

-G 

Sh 

ft) 

-G 

Ph 73 
ft) o 

g o 

G o 

u 

p 

ft) 

73 d _k> P 

2 ^ .2 G ft) 

^ ^ 73 q > 

00 ^ ftH ,P pQ 


VD 

O' 

oo 

o 

7 t 

CO 

7 t- 

in 

i—H 

O' 

7 j- 

oo 

04 

O' 

i—H 

OV 

o\ 

O' 

VD 

VD 

04 

O' 

O 

CO 

CO 

in 

i—H 

CO 

04 

04 

VD 

O' 

CO 

O 

05 

^ 05 

£ 

o 2 ° q 3 



O' 

OO 

VD 

04 

O' 

OV 

VD 

OV 

OO 

O' 

o 

in 

1 —H 

7 t 

VD 

VD 

OV 

oo 

04 

7 j- 

o 

i—H 

04 

O' 

m 

7 j- 

VD 

04 

m 

i—H 

VD 

o 

VD 

ftft 

ft) 

ft) V 

G 

o p o -9 P 


H 

7 t 

in 

o 

O' 

VD 

O' 

OO 

O' 


OV 

o 

04 

O' 

04 

O' 

o 

O' 

CO 

O 

0 ' 


O' 

CO 

in 

O' 

0 ' 

OV 


O' 

OO 

7 t- 

o 

o 

05 

G 

O 

ft) 

^ G 

O 

04 H -2 __i x 


o 

OV 

OV 

oo 

o 

O' 

o 

CO 

o 

oo 

i—H 

oo 

O' 

o 

O 

o 

oo 

oo 

7 |- 

OO 

o 

oo 

OO 

O' 

OV 

o 

in 

m 

O' 

o 

7 j- 

VD 

in 

O' 

05 


s « ^ bs 


P 4 

q 

OV 

OO 

OV 

o 

OV 

o 

m 

o 

OV 

OV 

OV 

OO 

o 

7 f 

o 

OV 

OV 

OV 

OV 

o 

OV 

OV 

OV 

o 

o 

OV 

OV 

OV 

o 

OV 

OV 

OV 

o 



X 


CO 

VD 

CO 

o 

CO 

o 

m 

o 

CO 

VD 

CO 

VD 

o 

CO 

o 

CO 

CO 

VD 

CO 

o 

CO 

CO 

CO 

04 

o 

CO 

CO 

CO 

o 

CO 

CO 

CO 

04 

pH 

0 /) 

ftG ’1 X 3 2 


u 

04 

i—H 

04 

o 

04 

o 

o 

o 

04 

i—H 

04 

i—H 

o 

o 

o 

04 

04 

i—H 

04 

o 

04 

04 

04 

o 

o 

04 

04 

04 

o 

04 

04 

04 

o 

G 

-G 

q o 



04 

VD 

04 

wo 

04 

wo 

VD 

7 f 

04 

VD 

04 

VD 

in 

VD 

in 

04 

04 

VD 

04 

wo 

04 

04 

04 

VD 

7 f 

04 

04 

04 

tI" 

04 

04 

04 

VD 

^G 

G 

G 

ft) 

H 1 > H ^p 

^ O c/ M 


VD VD O 

m 04 in 

O h (N 
O O O 

o o o 
q q q 

G G G 

O O O 

§§§ 
cn oo cn 

u u u 


VD o~ 
OV CSV 
04 04 

o o 
o o 

q q 

G G 

O O 

§§ 
cn oo 

u u 


tj- in 
i —h 04 
CO CO 

o o 
o o 

q q 

G G 

o o 

§§ 
cn oo 
U U 


* 

OO OV 

in o- 

CO CO 

o o 
o o 

q q 

G G 

o o 

s § 

cn oo 

U U 


§§ 
cn oo 
U U 


* 

VD OO OV 
m O h 
7 f ^J- VD 

o o o 
o o o 

q q q 

G G G 

O O O 

cn cn cn 

u u u 


d- o 

in vd 

VD VD 
O O 
O O 

q q 

G G 

o o 

§§ 
00 00 
u u 


O' O' 
VD ^H 
VD O' 

o o 
o o 

q q 

G G 

O O 

s§ 
00 00 
U U 


d- — d- 

O H C4 
O' OO OO 

o o o 
o o o 

q q q 

G G G 

O O O 

§§§ 
cn cn cn 

u u u 


oo O' 
04 CO 
OV O 
O -H 

o o 
q q 

G G 

O O 

s § 
00 00 
U U 


OO Tt 

co m 

o o 

o o 
q q 

G G 

o o 

00 00 
u u 


H O ^ 
co d- d- 


o o o 
q q q 

G G G 

O O O 

§§§ 
cn cn cn 

u u u 


O VD OO 
VD OV O 
H (N CO 

o o o 
q q q 

G G G 

O O O 

§§§ 
cn cn cn 

u u u 


VD 04 
OO CO 
OV i-H 

VD ^|- 
O -H 

q q 

G G 

O O 

00 00 
U U 


OC- 
G ' 

S 

o «. 


•c 

<D 


ft) 

p b 

a) 

G 

X 

O o 

X 

W) 

_G 


G 

D ^ 

H 

o 

X 

05 

05 

G 

G 

5h X 

<23 

73 G 

b ^ 

2 

73 

G 

G 

73 

73 

G 

G 

G 

S G 

G 

O 

G 

G ft) 


ft-H 

q 

2 p 

2 

C/3 

ft) 

2 £ 

=3. 

73 



O 


•G 73 ’ 
o 

s Q 


2 c 5 


^ cn 

cn 


:5 


> 

- 5 i) 


73 52 
G -2 
c 3 O 


ftp w 

H o 

O o 

Ph O 
O O 

U K 


ftp o 
op ^ 

G 

o 

G 

ft) 


" C/3 

-d 

P 

G 

O cO 

c 2 

.22 2 


ft) oo 

h h o 

ft) 73 

ft) -G o 

P H -g 

Gh 


Article number, page 6 of 27 






P. T. McGinnis et al.: Probing the inner disks of AA Tau-like systems in NGC 2264 



Mon-000604 


(e) 


JO-E4SOOOD 


JD-2450000 





JD-245QO0D 


JD-24S000O 


JD-2450000 


Fig. 3: Examples of different light curve morphologies found in CTTS in NGC 2264. (a) and (b) show periodic light curves attributed 
to circumstellar extinction, (c) shows an accretion burst dominated light curve, (d) shows a light curve dominated by aperiodic 
extinction, (e) shows a spot-like light curve, and (f) shows a light curve whose morphology is very complex, probably as a result of 
a combination of processes. 


Mon-1308 



Mon-250 



Fig. 4:1-band and CoRoT photometry, shifted in magnitude to coincide, for the AA Tau-like stars Mon-1308 and Mon-250. The high 
cadence space-based observations allow the AA Tau-like flux dips to be easily identified, while with even high quality ground-based 
photometry it is difficult to obtain enough information for the variability to be quantified. 


3.2. Veiling and variability 

3.2.1. Correlation between photometric variability and veiling 

Veiling occurs as a result of a strong UV and optical excess pro¬ 
duced in the accretion shocks on the stellar surface of CTTS. It 
is present in spectra when the hot spot caused by these accre¬ 
tion shocks is visible to the observer. This UV excess results in 
an added blue continuum to the photosphere’s spectra, resulting 
in somewhat shallower absorption lines than those from a pure 
photosphere. For the CTTS in NGC 2264 for which we obtained 
FLAMES spectra, we measured the veiling and its variability, us¬ 
ing the equivalent widths of the Li I 6707. 8A and Fe I 6633.4A 
lines. These lines were chosen since they are the absorption lines 


with the best signal-to-noise ratio in our spectra. As veiling in¬ 
creases, the equivalent widths we measure will decrease. 

Many of these CTTS show a certain amount of variability in 
their measured veiling. We verified whether there was a correla¬ 
tion between this variability and the photometric variability. In 
principle, a correlation or anti-correlation should exist in some 
cases, depending on the physical processes occurring on the stel¬ 
lar photosphere. If, for instance, the variability is due to obscura¬ 
tion of the photosphere by an inner disk warp that is associated 
with accretion columns, we can expect two possible scenarios. 
When the accretion shocks on the stellar surface are not visi¬ 
ble, or the veiling is too small to be accurately measured, there 
will be no considerable variability in the observed line equiva- 
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Fig. 5: Top: plots of LilEW vs. I magnitude for three different stars. Bottom: I-band and CoRoT light curves, shifted in magnitude 
to coincide, of the same three stars. Left panels: the star has an accretion-burst dominated light curve, and veiling variability that 
supports the variable accretion scenario. Middle panels: the star has a light curve that presents aperiodic extinction events, and 
veiling variability that supports a scenario in which the occulting structures are associated with accretion shocks. Right panels: the 
star has a light curve with A A Tau-like characteristics, while its near lack of veiling variability can neither confirm nor refute the 
hypothesis that it is an AA Tau-like system. 


lent widths. If, on the other hand, accretion shocks appear in our 
line of sight and produce measurable veiling, they should appear 
along with the inner disk warp, since both are associated with 
the accretion columns. In this case, when the inner disk warp 
occults the star causing brightness minima, the veiling should 
increase thanks to the appearance of the hot spots caused by ac¬ 
cretion shocks. Therefore, for an A A Tau-like light curve, the 
equivalent width of Li I (LilEW) or Fe I (FelEW) should either 
remain relatively constant or increase with increasing brightness. 
In the latter case, this could be evidence that the star’s photomet¬ 
ric variability is due to obscuration by an inner disk warp that is 
spatially associated with the accretion columns. 

For photometric variability that is due mainly to a config¬ 
uration of hot spots on the stellar surface, the veiling should 
increase with increasing brightness, and therefore FelEW and 
LilEW should decrease with increasing brightness. This is the 
case for the accretion-driven flux bursts discussed by Stauffer 
et al. (2014). The reoccurring narrow increases in brightness are 
likely due to an increase in accretion luminosity and are there¬ 
fore associated with an increase in veiling. 

If the star’s photometric variability is due simply to the rota¬ 
tional modulation of a stable, cold spot on the photosphere, the 
veiling would not be expected to vary, since cold spots do not 
cause veiling. 

Measures of veiling using the Fel 6633 A A line are much less 
accurate than Lil 6707.8A, since it is shallower in these stars and 
therefore suffers more from noise in the spectra. For this reason 
veiling was only measured using Fel for a small subset of these 
stars with FLAMES spectroscopy. In these cases the observed 
tendencies agree with those determined using the Lil line. 

We made plots of LilEW vs. USNO I-band magnitudes for 
the stars in our sample where both data were available. The I- 
band light curves were used rather than the CoRoT light curves 
because they had better time coverage, despite the lower ca¬ 


dence. There are at most six FLAMES spectra of each star ob¬ 
tained simultaneously with the CoRoT observations, whereas the 
USNO observations spanned the entire FLAMES campaign. In 
order to construct the plots, I-band measurements taken on the 
same night were averaged and matched with the equivalent width 
found using the FLAMES spectrum taken on that night. The 
USNO and VLT are separated by little more than 40° in lon¬ 
gitude, so the observations were likely taken within a few hours 
of each other on most nights. 

Fig. 5 shows plots of LilEW vs. I-band magnitudes (top pan¬ 
els) for three different examples, along with the corresponding 
CoRoT and I-band light curves (bottom panels). The left pan¬ 
els of Fig. 5 show a star whose veiling increases with increasing 
brightness. We can see from its light curve that it is dominated 
by flux bursts attributed to variable accretion events. The mid¬ 
dle panels of Fig. 5 show a star whose veiling increases with 
decreasing brightness, a tendency we attribute to occultation of 
the photosphere by an inner disk warp associated with the accre¬ 
tion mechanism. Its light curve morphology shows variability 
that could be explained by this scenario. The right panels of Fig. 
5 show a star whose veiling does not vary considerably, which 
can mean the photometric variability is due to a configuration of 
cold spots on the surface, or to obscuration by an inner disk warp 
that either fully occults the accretion shocks associated with it, 
leaving the veiling unaltered during eclipses, or the accretion is 
relatively low, and only a small hot spot is present, not capable of 
generating a large UV excess or veiling. The light curve agrees 
with the inner disk warp scenario, since cold spots would result 
in a much more regular photometric variability, but the veiling 
can neither confirm nor disprove this hypothesis. This star also 
presents little UV excess, consistent with either of these scenar¬ 
ios. 

The star Mon-811 was outside the USNO field, and therefore 
we have no I-band information for it. For this one star, we made a 


Article number, page 8 of 27 





























P. T. McGinnis et al.: Probing the inner disks of AA Tau-like systems in NGC 2264 


Mon-441 


Mon-498 



I Magnitude 
Mon-660 



I Magnitude 
Mon-1054 


0.50 

0.48 
' 0.46 
0.44 

0.42 



14.60 14.29 13.97 13.65 13.34 13.02 
I Magnitude 



I Magnitude 


Fig. 6: Plots of LilEW vs. I magnitude for some of the stars in our 
sample. The star Mon-1054 is an example of a candidate AA Tau 
system that shows increasing LilEW with decreasing magnitude. 
The star Mon-441 shows similar behavior, though its photomet¬ 
ric behavior is aperiodic. The star Mon-660 is an AA Tau-like 
system that appears to show increasing LilEW with decreasing 
magnitude, though the uncertainties are quite large, and the star 
Mon-498 is an example of a system with very irregular veiling 
that shows no obvious pattern. 


plot of LilEW vs. CoRoT flux counts instead. This plot was not 
nearly as well sampled, since very few FLAMES spectra were 
taken during the CoRoT observations, but it shows a tendency 
similar to the middle panels of Fig. 5. 

Fig. 6 shows plots of Lil equivalent width vs. I magnitude for 
some of the A A Tau-like or aperiodic extinction stars in our sam¬ 
ple (plots of the missing stars can be found in the online material, 
Fig. 21). In many cases the veiling variability is very irregular, 
and we cannot see a clear pattern (e.g., Mon-498). For the star 
Mon-379 (shown in the right panels of Fig. 5), the veiling is rel¬ 
atively constant. For six stars, we see an apparent decrease in 
veiling (increase in Lil equivalent width) as the star’s brightness 
increases, which could indicate that the mechanism responsible 
for their photometric variability is obscuration by an inner disk 
structure associated with accretion columns. This is the case for 
Mon-250, Mon-441, Mon-660, Mon-811, Mon-1054, and Mon- 
1144. The last has the largest variability amplitude in veiling 
(~ 0.3A), meaning there should be significant hot spots on the 
stellar surface. It is the star with the highest UV excess among 
the stars in our sample where veiling was measured, which could 
be expected, since the hot spots associated with the accretion 
shocks on the surface of the star should produce a considerable 
UV excess when visible. 

Two of the stars that show an increase in veiling during the 
minima in their light curves have aperiodic photometric variabil¬ 
ity (Mon-441 and Mon-1144). This provides support for our hy¬ 
pothesis that the occultations in these light curves are caused by 
dust lifted above the disk plane at the base of unstable accretion 
streams, since these should be associated with accretion shocks 
on the stellar surface in a manner analogous to the stable ac¬ 
cretion funnels. If these stars are undergoing unstable accretion, 
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Fig. 7: Mon-250, a) 2011 Corot light curve folded in phase, b) 
H a line profile during different rotational phases, indicated in 
each panel. A redshifted absorption component seen only near 
phase 0.5 is indicated by a red arrow, c) Measured veiling folded 
in phase, d) Shift in radial velocity measured from photospheric 
lines in the spectra, folded in phase. In all four panels, the same 
period of 8.6 days, determined from the radial velocity variation 
of photospheric lines, was used to calculate the phase. The same 
initial date was also used. 


they should show no obvious periodicity in their light curves, but 
if they are seen at high enough inclinations, they should share 
other similarities with A A Tau-like systems, such as the one we 
see here. 

3.2.2. The specific case of Mon-250 

We analyzed the H a line profile of the star Mon-250 during dif¬ 
ferent rotational phases (Fig. 7b), and note a redshifted absorp¬ 
tion component (indicated by an arrow in the figure) that ap¬ 
pears in several spectra during the light curve minima that is not 
present in other phases. This is clear evidence that the accretion 
funnel and hot spot appear in front of our line of sight during the 
photometric flux dips, since the infalling material responsible for 
the absorption is moving away from us. 
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This star’s veiling variability also points to the appearance 
of accretion shocks during occultations. For this particular star, 
we were able to measure veiling in all 20 FLAMES epochs by 
comparing line depths to those of a weak line T Tauri star of the 
same spectral type, to determine what fraction of the observed 
continuum should be due to the hot spot. Fig. 7c shows a plot 
of veiling folded in phase, using the rotation period determined 
from radial velocity shifts of photospheric lines in the spectra 
(Fig. 7d). It can be seen that there is almost always an increase 
in veiling as the light curve undergoes a minimum in brightness. 
This is evidence of the appearance of a hot spot on the stellar 
surface, one large enough to produce measurable veiling, dur¬ 
ing occultations. The rotation period found from the variability 
of radial velocities measured in the Mon-250 photospheric lines 
was of 8.6 ± 0.5 days, consistent with the period found from the 
CoRoT light curve observed at the same epoch of 8.93 ± 0.50 
days. This shows that the occulting structure is located at or near 
the disk’s co-rotation radius, as was the case for AA Tau. Un¬ 
fortunately, this type of analysis was only possible for this star 
since it was the only one with a high enough variability ampli¬ 
tude in veiling and radial velocity of photospheric lines to draw 
significant conclusions. 

We attribute these phenomena, which were also observed for 
A A Tau (Bouvier et al. 2007), to the appearance of a major accre¬ 
tion funnel and accretion shock associated with it, at the moment 
of occultation by the inner disk warp. This shows the association 
between the inner disk warp and the stable accretion columns, 
and provides support to our hypothesis that in these cases, the 
same interaction between the inclined stellar magnetic field and 
the inner disk that is causing the material to accrete onto the star 
via two major accretion columns is also responsible for lifting 
dust from the disk plane, creating a warp in the disk near the 
co-rotation radius. 


3.3. Stable and unstable accretion regimes 

In this section, we discuss the main differences between the 
quasi-periodic AA Tau-like stars and the stars whose light curves 
are dominated by circumstellar extinction, but have no obvious 
periodicity. We propose that the former may be undergoing ac¬ 
cretion via a stable regime, while the latter could be accreting 
via an unstable regime. Kurosawa & Romanova (2013) cite the 
parameters that are important in determining through which of 
these regimes a star will accrete. They are mass accretion rate 
(M acc ), truncation radius, stellar rotation rate, the inclination an¬ 
gle /3 between the stellar rotation axis and the axis of the mag¬ 
netic field (Q and p, respectively, in Fig. 1), and the tr-parameter, 
which describes viscosity in the disk. 

Kurosawa & Romanova (2013) mention a possible transition 
between regimes, which can occur if one or a few of the factors 
that determine instability change with time within a certain sys¬ 
tem. For instance, if the truncation radius R mag becomes smaller 
than the co-rotation radius R co , due for example to a change in 
the stellar magnetic field, accretion flows tend to become unsta¬ 
ble. The mass accretion rate can also affect the ratio R mag /R co , 
since a higher mass accretion rate leads to a smaller trunca¬ 
tion radius, therefore leading to a more unstable configuration. 
The viscosity parameter (a) controls the global mass accretion 
rate. The higher a is, the higher the mass accretion rate is, lead¬ 
ing to more instabilities. According to Kulkarni & Romanova 
(2009), the transition from a stable to an unstable regime occurs 
around critical values of mass accretion rates. They state typi¬ 
cal critical values for stars with R mag « 2 - 3 R* to be around 


M acc = 1.7 x 10 7 , and for stars with R mag « 4 - 5R*, around 
M acc = 2.1 x 1(T 8 . 

A smaller angle p also leads to less stable configurations, 
since at large ft the magnetic poles become closer to the inner 
disk rim, lowering the potential barrier in the vertical direction 
that the material must overcome to be lifted into the accretion 
funnels. According to Kulkarni & Romanova (2009), values of 
P > 25° usually lead to stable configurations (see also Romanova 
et al. 2014). 

Venuti et al. (2014) showed that mass accretion rates in 
CTTS can vary considerably over time. They calculated mass 
accretion rates for stars in NGC 2264 using data from 2010 De¬ 
cember and 2012 January, the latter being nearly simultaneous 
with the CoRoT observations. They show that many of the stars 
in our sample had mass accretion rates vary by up to 50% in little 
over one year. Unfortunately, there is no estimate of mass accre¬ 
tion rate simultaneous with the 2008 CoRoT observations, but 
it is possible that a variation such as this contributed to a tran¬ 
sition between accretion regimes. However, Venuti et al. (2014) 
also show that long term changes in mass accretion rates are not 
necessarily larger than short term changes. 

Magnetic field configurations may also be responsible for 
these transitions. It is possible for the magnetic field of a CTTS 
to change significantly in a few years. This has been observed, 
for example, in the CTTS V2129 Oph (Donati et al. 2011). Its 
magnetic field was measured in two different epochs, nearly four 
years apart, and the magnetic field strength varied greatly from 
one epoch to the other. The dipolar component was found to 
be about 3 times stronger in the second epoch, and the octupo- 
lar component 1.5 times stronger than four years earlier. Donati 
et al. (201 1) found that the ratio R mag /Rco nearly doubled in that 
time. 

In this study we have seen that 38%+\\% (5/13) of the stars 
that were observed in 2008 as having photometric variability 
similar to AA Tau were observed again in 2011, showing ape¬ 
riodic variability due to extinction from circumstellar material. 
Of the stars that were observed in 2008 as having aperiodic 
photometric variability due to circumstellar extinction, 62%+}^ 
(8/13) were observed again in 2011, then showing variabil¬ 
ity similar to AA Tau. If our proposed scenario is true, then 
this points to a transition from a stable to an unstable accre¬ 
tion regime, or vice-versa, in a matter of less than four years, 
for 51% of our sample of possible AA Tau-like stars, i.e., 
(5+8)/23. An example of two of the stars that appear to have 
undergone this transition is shown in Fig. 8. 

We also performed a periodogram analysis of the light curves 
of these 13 stars that changed between periodic and aperiodic 
photometric behavior. The periodograms of each epoch were an¬ 
alyzed separately and compared. For 11 stars, during the epoch 
of unstable accretion at least a weak signal was present at the 
period found for the AA Tau phase, believed to be the stellar 
rotation period. For eight of these, this peak was the only one 
present, or the strongest, despite being much weaker than the 
one present during the A A Tau phase. One exception is the star 
Mon-441, which shows a weak signal during its unstable phase, 
but this signal does not correspond to the period found in the AA 
Tau phase, or to a multiple of it. 

In most cases the variable depths of the flux dip amplitude 
in a star’s light curve during its aperiodic phase are smaller than 
those measured in its AA Tau phase by 20% to 60%. However, 
for 4 of the 13 stars that underwent a transition between AA 
Tau-like and aperiodic light curve, the flux dip depths remained 
nearly the same in both epochs, and for 2 more stars the variabil¬ 
ity amplitude during their aperiodic phase was larger by 20% to 


Article number, page 10 of 27 


P. T. McGinnis et al.: Probing the inner disks of AA Tau-like systems in NGC 2264 


Mon-297 Mon-297 



JD-2450000 JD-2450000 

Mon-774 Mon-774 




JD-2450000 JD-2450000 

Fig. 8: Corot light curves of Mon-297 and Mon-774 during their 
AA Tau phases (left) and aperiodic phases (right). 


30%. This means that, at least in some cases, during the unstable 
accretion phase dust must be lifted above the disk plane at sim¬ 
ilar, and sometimes even higher, altitudes as during the stable 
accretion phase. 

As for the star Mon-928, the flux dips present in 2008 disap¬ 
pear entirely in 2011. The light curve we observe in the second 
epoch is not characteristic of circumstellar extinction. The pe¬ 
riod of 9.92 days found during the first epoch is very similar to 
periods previously measured for this star in the literature (Lamm 
et al. 2005), yet it is non-existent in the 2011 light curves (see 
Fig. 9). One possibility is that in 2011 the inner disk warp did not 
reach high enough altitudes above the disk to occult the stellar 
photosphere. At the inclination of 71° found by us for this sys¬ 
tem (see Sect. 3.5), our simulations using the occultation model 
show that if the warp found to best reproduce the 2008 CoRoT 
light curve decreases in height by at least 26%, the occultations 
observed in 2008 would no longer occur. This amount of varia¬ 
tion in warp height is well within the range we observe for these 
objects (see Sect. 3.5). This star is also one of the two that appear 
to have undergone a small increase in mass accretion rate from 
2010 December to 2012 January (Venuti et al. 2014); therefore, 
it is possible that a higher mass accretion rate led to a more un¬ 
stable configuration during the 2011 observations, resulting in a 
less pronounced disk warp. 

Kurosawa & Romanova (2013) also predict that for a more 
unstable regime, one should observe a higher frequency of flux 
peaks per stellar rotation period in the light curves of stars of 
medium or low inclinations thanks to the presence of a larger 
number of accretion shocks on the photosphere. It is fair to as¬ 
sume that, for systems seen at high inclinations, a higher fre¬ 
quency of flux dips would also be observed for more unsta¬ 
ble scenarios because of the larger number of accretion streams 
present around the star, whose base could contain dust that can 
occult the photosphere. Of the stars in our sample that changed 
from one regime to the other, most show two to three flux dips 
per stellar rotation period during the unstable accretion phase, 
while some show up to seven flux dips per rotation period. This 
may be an indication of the degree of instability of each star-disk 


Mon-928 Mon-928 



Fig. 9: CoRoT light curves of Mon-928 in 2008 (left) and 2011 
(right). 

system. The two stars that show the highest frequencies of flux 
dips during their unstable accretion phases are the two that have 
the lowest masses in our sample of possible AA Tau-like stars. 
As is discussed in Sect. 3.7, low mass stars appear to present 
more unstable accretion, possibly as a result of magnetic field 
configurations. These higher frequencies of flux dips can also be 
caused by dust lifted above the disk plane because of magne- 
torotational instabilities that do not necessarily lead to accretion 
streams (Turner et al. 2010; Hirose & Turner 2011). 

Of the 23 stars that present AA Tau-like behavior, only six 
maintain this behavior in both 2008 and 2011. Many appear to 
show some signs of weak instability even during their AA Tau 
phase, sometimes presenting two or three flux dips in one rota¬ 
tion cycle. A few examples are Mon-297, Mon-456, and Mon- 
824 (Fig. 10). Most show many small secondary flux dips present 
in the light curve, indicating a possible coexistence of secondary 
accretion streams with the main accretion funnels. 

3.4. Comparing optical and infrared light curves 

Of the 33 stars in our sample of A A Tau-like or aperiodic extinc¬ 
tion dominated CoRoT light curves, 29 were observed almost si¬ 
multaneously with Spitzer in mapping mode. Fig. 11 shows a few 
examples of Spitzer IRAC 3.6/mi and 4.5yum light curves, plot¬ 
ted along with the respective CoRoT light curves. The CoRoT 
flux was transformed into magnitude using its greatest value as 
magnitude zero, then shifted to overlap the Spitzer IRAC data. 
A figure showing the light curves of all other stars in our sample 
can be found in the online material, Fig. 22. 

We see that 4 of the 29 stars show very similar behavior in 
both bands, and present A A Tau-like modulation in the infrared 
as well as the optical. These stars are Mon-660, Mon-811, Mon- 
1140, and Mon-1308. Another five show similar behavior, with 
minima and maxima generally coinciding in the infrared and op¬ 
tical, but these light curves do not overlap perfectly. These are 
Mon-297, Mon-314, Mon-433, Mon-441, and Mon-456. 

There are nine stars that correlate in some parts of their in¬ 
frared and optical light curves, while other parts behave dif¬ 
ferently, or even anti-correlate at times. Stars Mon-56, Mon- 
126, Mon-358, Mon-654, Mon-717, Mon-774, Mon-1054, Mon- 
1144, and Mon-1167 show this behavior. The 10 stars Mon-296, 
Mon-325, Mon-379, Mon-498, Mon-619, Mon-667, Mon-928, 
Mon-1037, Mon-1038, and Mon-1131 have entirely different in¬ 
frared and optical light curves, with very little or no features in 
common. 

The star Mon-250 has Spitzer IRAC light curves that almost 
anti-correlate with its CoRoT light curve. Most of the Spitzer 
IRAC maxima coincide with CoRoT minima, meaning this vari- 
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Fig. 10: CoRoT light curves of Mon-297, Mon-456, and Mon-824 during their AA Tau phase. These stars sometimes present two, or 
even three, flux dips per rotation period, and often secondary flux dips are present alongside the main flux dips, indicating a possible 
coexistence between unstable accretion streams and stable accretion funnels. 
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Fig. 11: Spitzer IRAC 3 .6pm and 4.5pm light curves (blue and red filled circles, respectively) overplotted on CoRoT light curves 
(black circles). The CoRoT and IRAC 3 .6pm light curves were shifted in magnitude for easier comparison. 


ability cannot possibly be the consequence of a hot or cold spot 
on the stellar surface, since any variability on the surface should 
be present in ah wavelengths. To explain why there is more IR 
emission during the optical eclipses, we propose the following 
scenario. The stable accretion regime assumes that there are two 
main accretion funnels opposite each other, one in each hemi¬ 
sphere. Each of these funnels is associated with a hot spot on the 
stellar surface, where the accretion shock is located. The base 
of each accretion funnel is associated with part of the inner disk 
warp. When one part of the warp occults the stellar photosphere 
in our line of sight, the hot spot opposite it, as well as the rest 
of the stellar photosphere on that side of the star, illuminates the 
corresponding part of the inner disk warp. The dust located there 
absorbs the light from the hot spot and photosphere and re-emits 
it at longer wavelengths, causing an increase in emission in the 
infrared (Fig. 12). This emission in the infrared is predicted by 
the radiation transfer model of Whitney et al. (2013) when ap¬ 


plied to an accretion hot spot associated with an inner disk warp. 
This behavior is also observed in a few other light curves, dur¬ 
ing one or two rotation cycles (see, e.g., Mon-1054 in Fig. 11, 
between MJD 55909 and MJD 55913). 

The aiRAc index is a good indicator of the amount of dust 
present in the inner accretion disk. It corresponds to the slope 
of the spectral energy distribution (SED) between 3.6/im and 
8/im. Unfortunately, there were no Spitzer IRAC observations 
in 8pm simultaneously with the 2011 CoRoT-Spitzer campaign, 
since by this time Spitzer was already operating in its warm mis¬ 
sion. We do, however, have Chirac values from Teixeira et al. 
(2012) for 14 of these 29 stars with simultaneous CoRoT and 
Spitzer light curves, and archival IRAC and Multiband Imaging 
Photometry for Spitzer (MIPS) data (for details see Cody et al. 
2014), which enabled us to calculate values of Chirac for the stars 
that were not included in Teixeira et al. (2012). These values are 
given in Table 1 . Lada et al. (2006) established a criterion to clas- 
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Fig. 12: When one side of the inner disk warp occults the stellar 
photosphere, the hot spot and photosphere on the opposite side 
illuminate the part of the warp on that side. The light is absorbed 
by the dust and is re-emitted at longer wavelengths. 


sify the inner disk structures of systems according to the Chirac 
index. Stars with a IR Ac < -2.56 are classified as naked photo¬ 
spheres (systems devoid of dust within a few 0.1 AU), those with 
-2.56 < Chirac < -1.80 have anemic disks (their inner disk re¬ 
gion is optically thin), those with -1.8 < a IR Ac < -0.5 have 
optically thick inner disks, stars with -0.5 < Chirac < 0.5 are flat 
spectra sources, and stars with Chirac > 0-5 are Class I protostel- 
lar objects. 

We note that the objects in our sample whose CoRoT and 
Spitzer light curves correlate best have values of Chirac between 
-1.4 and -1.9, meaning they are generally in the transition region 
between optically thin and optically thick inner disks. If there 
is too much dust within a few AU of the star (higher values of 
Chirac )»the Spitzer light curves would be dominated by emission 
from warm circumstellar dust, rather than eclipses of the stellar 
photosphere. On the other hand, if the portion of the disk oc¬ 
culting the star is not optically thick, then the flux dips at IRAC 
wavelengths could be much shallower than for the case where 
we assume optically thick absorption at all bands. This seems to 
be the case for Mon-296, whose Spitzer light curve remains rela¬ 
tively steady except for a very shallow dip occurring at the same 
time as the optical light curve’s lowest minimum, and a small 
burst near one of its maxima (top right panel in Fig. 11). This 
star’s Chirac shows that it has an anemic disk. The star Mon-250, 
mentioned in the discussion above, also has an anemic disk. It 
has a value of Chirac = -2.10, one of the lowest in our sample. 

Another useful quantity is the ratio of the disk flux to the 
stellar flux at 4.5 pm, estimated from the slope of each star’s 
SED as compared to an estimate of the photospheric flux based 
on the star’s spectral type and the Kurucz-Lejeune models (Leje- 
une et al. 1997), normalized to the observed J -band and assum¬ 
ing an appropriate amount of reddening customized to each ob¬ 
ject. This indicates what fraction of the 4.5 pm flux originates at 
the disk, and should relate directly to how well the light curves 
correlate. If this ratio is zero, then the optical and infrared light 
curves should coincide. If it is considerably larger, very little of 
the infrared flux comes from the star, and the flux dips present in 
the optical light curves should be absent from the infrared light 
curves. The values of F 4< 5 yUm (disk)/F 4< 5 yUm (star) estimated for our 
stars are shown in Table 1 . 

We see some of these predicted tendencies, but the relation¬ 
ship is not as direct as one would expect. The stars with the 
lowest fractions of disk flux at 4.5 pm tend to have well cor¬ 
related light curves, and generally very similar variability ampli¬ 
tudes in the optical and infrared, while those with a very high 
value of disk flux per stellar flux have uncorrelated, or weakly 
correlated, optical and infrared light curves. For the intermedi¬ 
ate values, however, we find all sorts of scenarios. Mon-660, 
for instance, has an estimated ratio of disk flux to stellar flux 


of F 4 . 5 /im (disk)/F 4 . 5 /im (star) = 3.4, and very strongly correlated 
light curves, while for Mon-928 F 4 . 5 yUm (disk)/F 4 . 5 yUm (star) = 1.4, 
and its optical and infrared light curves are completely different. 
This could be due to uncertainties in the SED fits caused by poor 
or time variable data, or could simply be due to the time vari¬ 
able nature of the SEDs. There could have been some changes in 
the disk emission between the Spitzer cryogenic mission, from 
whence the values used to plot the SEDs were taken, and the 
2011 CoRoT and Spitzer observing campaign. 

For the stars where a correlation can be identified, we calcu¬ 
lated which extinction law is necessary to transform the optical 
to the infrared light curve, assuming an optically thin regime. Ex¬ 
tinction laws Ax/Ar , where A = 3.6pm and 4.5yum, were found to 
be between 0.28 and 0.77. According to Cardelli et al. (1989), we 
can derive values of Ax I Ay for a given R v (where R v is defined as 
A V /E(B- V )). A value of R v = 3.1, typically used for interstellar 
matter (ISM), translates to A^/Ay = 0.05 and A^.s/Ay = 0.04. 
Transforming to A R , since the CoRoT filter is centered close to 
the center of the R Johnson filter, we would have A^/Ar = 0.07 
and A 4< 5 /A r = 0.05. These extinction laws are much lower than 
those found for the stars in our sample, implying that the dis¬ 
tribution of dust grains in the inner circumstellar disk region is 
quite different from ISM, possibly containing fewer small grains. 
In addition, the observed dip depth ratio should be corrected for 
the fraction of IR light emitted by the disk. This would only in¬ 
crease our values of Ax/A R . 

3.5. Modeling the CoRoT light curves 

In this section we use the occultation model proposed for AA 
Tau in Bouvier et al. (1999) (described in Sect. 1.1) to attempt 
to reproduce the A A Tau-like light curves in our sample. If the 
optical photometric variability observed is due to periodic occul- 
tations by an inner disk warp, as we propose, then this geomet¬ 
rical model should be capable of reproducing the general trend 
of these light curves. We vary the model’s main parameters in an 
attempt to obtain the configuration that most closely matches the 
observed variability. In this way we determine which values of 
warp height and azimuthal extension are needed to explain these 
stars’ variability and how they must change from one rotation 
cycle to another in order to account for the changes in width and 
amplitude of each flux dip (see also Fonseca et al. 2014, for a 
detailed study of the light curve of V354 Mon - Mon-660). 

As was mentioned in Sect. 1.1 , the parameters of the occul¬ 
tation model are the star-disk system’s inclination i, the radius 
R w at which the warp is located, and the warp’s maximum height 
h max and azimuthal extension cp w . 

We assume that the periods we measured from the CoRoT 
light curves are the Keplerian periods of the occulting material 
orbiting the star; therefore, if we know the star’s mass, we can 
use Kepler’s third law in the following form to calculate the ra¬ 
dius at which the warp is located, 

(2) 

where M* is the star’s mass, Pk ep is the Keplerian rotation period, 
and G is the gravitational constant. To calculate the inclination 
we use the following relation between vsin i, stellar radius R *, 
and stellar rotation period P rot : 

. . 2 nR„ . . 

vsmz = -sin i. (3) 

Trot 
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We used stellar masses and radii determined by Venuti et al. 
(2014), who compared effective temperatures and bolometric lu¬ 
minosities placed on the Hertzsprung-Russell diagram to pre- 
main sequence model grids of Siess et al. (2000). The values are 
shown in Table 2, along with our values of Keplerian periods 
found from CoRoT light curves in each epoch, and the corre¬ 
sponding Keplerian rotation radius. 

Values of the stars’ projected rotational velocities (vsin i) 
were obtained comparing FLAMES spectra, when available, and 
Hectochelle spectra otherwise, to synthetic spectra using the 
program Spectroscopy Made Easy (SME, Valenti & Piskunov 
1996) and Oleg Kochukhov’s Binmag3 2 , as well as the spec¬ 
trum synthesis code Synth3 (Kochukhov 2007) and atomic line 
files extracted from the Vienna Atomic Line Database (VALD) 3 
(Piskunov et al. 1995; Ryabchikova et al. 1997; Kupka et al. 
1999, 2000). We used photospheric lines to estimate the values 
of v sin i shown in Tables 1 and 2. For a few of our stars either no 
spectrum was available, or the signal-to-noise ratio of the pho¬ 
tospheric lines was insufficient for this analysis and therefore no 
value is given. 

Co-rotation radius A stellar rotation period independent of 
CoRoT light curves was measured for only one star in our sam¬ 
ple, Mon-250, by determining the variability of the radial veloc¬ 
ity of photospheric lines in the FLAMES spectra (Fig. 7). This 
method was unsuccessful for the other stars, since their radial ve¬ 
locities did not present significant variability within our errors. 
The spectroscopic period found for Mon-250 was P rot = 8.6+0.5 
days, a value consistent with the periods found from its CoRoT 
light curves of P^ ep = 8.3 ± 1.0 days in 2008 and P^ ep = 8.9 ± 0.5 
days in 2011. The same period of 8.6 days was found by Sousa 
et al. (in prep.) in both the F la and the Hel 6678A lines. This 
indicates that the structure responsible for the occultations in the 
CoRoT light curves is located at or near the co-rotation radius. 

This was found to be true for the star AA Tau as well, when 
Bouvier et al. (2007) measured its rotational period using the ra¬ 
dial velocity of the Hel 5876A emission line and photospheric 
lines and found a value consistent with their photometric pe¬ 
riod. Therefore we find it reasonable to assume that the inner 
disk warp is located at the disk’s co-rotation radius for the other 
A A Tau-like systems. Since we do not have an accurate mea¬ 
sure of stellar rotation periods for the stars in our sample besides 
Mon-250, we use the periods determined from their CoRoT light 
curves in Equation 2 to calculate their inclinations. 

In order to further investigate whether these Keplerian pe¬ 
riods can also represent the stars’ rotation periods, we made a 
period histogram of all of the CTTS in NGC 2264 with CoRoT 
light curves classified as spot-like or AA Tau-like (see Fig. 13). 
The spot-like light curves are due to rotational modulation of sta¬ 
ble configurations of spots on the stellar surface, and therefore 
these values represent the stars’ rotational periods. The values 
from AA Tau-like light curves represent the Keplerian rotation 
periods of the inner disk material that occults the stellar photo¬ 
sphere. We see that, though the distributions are not identical, the 
values for AA Tau-like light curves fall within the region where 
rotational periods are found. 

Inclination We calculate sin i of each star using Equation 3 and 
taking the period found from the CoRoT light curve as P rot , ex¬ 
cept in the case of Mon-250, for which the spectroscopic period 

2 http://www.astro.uu.se/~oleg/ 

3 http://vald.astro.uu.se/ 



period (days) 

Fig. 13: Period histogram of classical T Tauri stars in NGC 2264 
whose CoRoT light curves were classified as AA Tau-like (blue) 
or spot-like (red). The periods of spot-like light curves represent 
stellar rotational periods, while those of A A Tau-like light curves 
represent Keplerian rotation periods at the warp location. The 
periods of spot-like light curves were taken from the 2011 data, 
while AA Tau-like light curves from both 2008 and 2011 were 
used. 


was used. For those stars observed to be periodic in both epochs, 
but for which a slightly different period was found in 2008 and 
2011, the value with the smallest error was used, i.e., the one 
whose periodogram presented a peak with smaller width at half 
maximum and higher intensity. The results are shown in Table 2. 

It is important to note that the inclination is a very difficult 
parameter to determine precisely, especially when dealing with 
high inclinations. In these cases a relatively small error in sin i 
is translated into a large error in i. Unfortunately our errors in 
sin i themselves are relatively large to begin with, therefore our 
values for i are very imprecise. 

Hartmann (2001) and Rebull et al. (2002) discuss the dif¬ 
ficulties in obtaining precise values of luminosity, temperature, 
and radius for TTSs because of their strong variabilities, spec¬ 
trum veiling, surface spots, and uncertainties in calculating the 
distance. The former showed that typical errors in radius deter¬ 
mination are on the order of 20%, while the latter show errors 
of typically 12% for TTSs in the spectral range K5-M2. In this 
study we consider errors of 20% in radius. 

Errors in period and v sin i vary from one star to another. For 
periods, we adopted the half width at half maximum of the peak 
in each periodogram as the uncertainty, which is typically around 
5% to 7%. Our errors in v sin i are typically on the order of 10% 
- 12%. This leads to errors around 35% in sin i, which is rather 
high, but consistent with other studies of TTS inclinations (e.g., 
Artemenko et al. 2012). 

Although the actual values of inclination are imprecise, it is 
not difficult to distinguish between systems with low inclinations 
and those with medium to high inclinations, since for low incli¬ 
nations errors in sin i do not translate into such large errors in i 
as with high inclinations. We see from Table 2 that the star Mon- 
1131 is seen at a low inclination, and we can consequently rule 
out the proposed occultation mechanism as the primary cause for 
the light curve variations. However, it is possible that there are 
two structures in each hemisphere, separated by 180°, obscuring 
the photosphere and making it seem as if the period were half the 
true value. If this were the case its inclination would be 64° 
though further investigation is needed to determine its true rota¬ 
tion period. For the time being this star has been excluded from 
further investigation using the occultation model, though it is 
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Table 2: Stellar parameters initially calculated for the A A Tau-like candidates. 


CSIMon ID 

Lbol 

(L @ ) 

M* 

(M 0 ) 

R, 

(Re) 

p i 

(days) 

Pi 

(days) 

v sin i 
(km/s) 

sin i 

Inclination 

o 

Rw, i 

0 R *) 

Rw,2 

(R„) 

CSIMon-000056 

0.77 

1.16 

1.52 

5.71 

5.86 

11.1+2.2 

0.85+0.38 

5+f 

9.29 

9.45 

CSIMon-000250 

1.23 

1.35 

1.63 

8.32 

8.93 

9.8+2.2 

1.02+0.49 

62±28 

11.71 

12.28 

CSIMon-000296 

1.57 

1.42 

1.71 

- 

3.91 

17.3+2.6 

0.78+0.33 

5ltH 

5.11 

6.86 

CSIMon-000297 

1.57 

1.42 

1.71 

3.16 

- 

30.8+1.8 

1.13+0.40 

69±21 

5.95 

- 

CSIMon-00035 8f 

0.22 

0.29 

1.39 

- 

5.86 

- 

- 

- 

- 

6.51 

CSIMon-000379 

2.12 

1.60 

1.99 

- 

3.68 

25.0+1.2 

0.91+0.30 

66!i 

5.95 

5.89 

CSIMon-000441 f 

0.50 

0.36 

1.92 

4.06 

- 

- 

- 


3.97 

- 

CSIMon-000456 

1.58 

1.41 

2.00 

- 

5.03 

19.3+1.5 

0.96+0.21 

74 + 33 

- 

6.92 

CSIMon-000498 

3.85 

1.90 

2.88 

4.23 

4.28 

24.1+2.4 

0.71+0.26 

4571 

4.73 

4.77 

CSIMon-000654 

0.47 

0.30 

2.03 

4.66 

- 

18.0+3.0 

0.82+0.37 

557| 

3.87 

- 

CSIMon-000660 

1.36 

1.40 

1.86 

5.25 

5.25 

21.4+2.5 

1.16+0.44 

68±22 

7.64 

7.50 

CSIMon-000774 

3.10 

1.83 

2.49 

3.46 

- 

30.7+1.3 

0.84+0.28 

57711 

4.73 

- 

CSIMon-000811 

1.08 

0.91 

1.97 

- 

7.88 

13.5+1.7 

1.07+0.35 

64+26 

- 

8.19 

CSIMon-000824 

1.96 

1.48 

2.23 

7.05 

- 

14.8+1.4 

0.92+0.40 

68 + 37 

7.90 

- 

CSIMon-000928 

0.58 

0.63 

1.65 

9.92 

- 

20.1+2.8 

2.38+1.05 

90 

10.09 

- 

CSIMon-001054 

0.63 

0.36 

2.17 

- 

4.08 

20.4+2.0 

0.76+0.28 

497g 

- 

3.52 

CSIMon-001131 

0.55 

0.36 

2.01 

- 

5.18 

8.8+1.7 

0.45+0.21 

27+14 

4.44 

4.46 

CSIMon-001140 

1.10 

1.31 

1.67 

3.87 

3.90 

19.8+2.3 

0.91+0.34 

667|f 

6.79 

6.83 

CSIMon-001167 

0.56 

0.30 

2.23 

- 

8.78 

11.3+2.1 

0.88+0.42 

627| 

5.18 

5.37 

CSIMon-001296 

0.91 

0.69 

1.99 

- 

9.75 

11.0+1.6 

1.07+0.48 

63 ±27 

- 

8.52 

CSIMon-001308 

0.53 

0.63 

1.59 

6.45 

6.68 

9.9+2.0 

0.82+0.37 

557| 

7.86 

8.04 

CSIMon-014132f 

0.37 

0.28 

1.88 

- 

4.48 

- 

- 


- 

3.98 


Notes. PI and P2 are the Keplerian periods found from the CoRoT light curves in 2008 and 2011, respectively, while R w , i and R w2 are the 
co-rotation radii found using PI and P2, respectively. The stars identified with a t either have no observed spectrum, or the S/N is too low to 
determine v sin /, and so we cannot calculate their inclinations. 


studied in Stauffer et al. (2015), where different mechanisms are 
proposed to explain its variability. 

For the other stars in our sample for which an inclination can 
be estimated, they tend toward high values, at least within the 
uncertainties (see Table 2). There are many stars for which the 
value of sin i found was larger than 1. In some cases a reasonable 
interval of values within the error bars falls under 1, in which 
case we argue that the system’s inclination should be anywhere 
between the inverse sine of the lowest value (sin i - A sin i, where 
A sin i is the uncertainty in sin i) and 90°. We then consider the 
inclination to be the mean value in this interval (for example, a 
star of sin i = 1.07 ± 0.32 has an inclination of 69° ± 21°, since 
the inverse sine of the lowest value, sin i = 1.07 - 0.32 = 0.75, 
is 48°). 

For the star Mon-928, its value of sin i is much larger than 
1, even when considering the error bars. It is possible that the 
material occulting this star is located farther away from the star 
than its co-rotation radius, in which case the Keplerian period we 
calculated using the CoRoT light curves is larger than the stel¬ 
lar rotation period, resulting in a larger value of sin i. It is also 
possible that the true period is half the value used. There is an ad¬ 
ditional, smaller, flux dip in its 2008 light curve located between 
the second and third larger flux dips (see Fig. 9), which justifies 
using a period of 4.96 days, rather than 9.92 days, though we 
gave preference to the latter since it coincides with the period 
found in the literature of 9.48 ± 2.34 (Lamm et al. 2005). If the 
period were indeed half the value we considered, this would lead 
to an inclination of 66° ± 24°. In either case, it is probably safe 
to consider that this system has a high inclination. Therefore, we 
used the occultation model described in Sect. 1.1 to find which 
interval of inclinations results in the best fit for the light curve, 


Mon-1308 


Mon-1308 



JD-2450000 JD-2450000 


Fig. 14: Fitted light curves of the star Mon-1308 using an in¬ 
clination of 55° (left) and 75° (right). In both cases the values 
of warp maximum height and azimuthal extension were set to 
those that best fit the width and amplitude of the minima without 
extrapolating to unphysical values. 


and consider this as the system’s inclination. We did this as well 
for the three stars for which we have no measure of v sin i, where 
we assume without proof that the systems are seen nearly edge- 
on. The new values are presented in Table 3. 

For the rest of the A A Tau-like stars, we used the occulta¬ 
tion model in a similar fashion as a fine tuning for the inclina¬ 
tion. We noticed that in many cases we were not able to suc¬ 
cessfully reproduce the observed widths in flux dips using the 
whole range of values of i initially calculated, unless we used 
other parameters that were unphysical (such as azimuthal ex¬ 
tension > 360°). Therefore, we took all values of i within 
the error bars and chose only the interval that could successfully 
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Table 3: Occultation model parameters found for the A A Tau-like candidates. 


CSIMon ID 

Epoch 

R W (R*) 

hmax/Rw 


Inclination (°) 

CSIMon-000056 

2008 

9.29 

0.24-0.25 (±0.10) 

320 (±40) - 360 

71±6 

CSIMon-000056 

2011 

9.45 

0.24-0.26 (±0.10) 

290 (±70) - 360 

71±6 

CSIMon-000250 

2008 

11.71 

0.23-0.28 (±0.05) 

140 (±30) - 360 

74±3 

CSIMon-000250 

2011 

12.28 

0.24-0.27 (±0.05) 

270 (±30) - 360 

74±3 

CSIMon-000296 

2011 

6.86 

0.21-0.31 (±0.09) 

290 (±40) - 360 

72±5 

CSIMon-000297 

2008 

5.95 

0.17-0.28 (±0.04) 

180 (±30)- 360 

75±2 

CSIMon-000358 

2011 

6.51 

0.17-0.19 (±0.06) 

210 (±50) - 360 

73±4 

CSIMon-000379 

2011 

5.89 

0.22-0.34 (±0.12) 

110 (±30)-270 (±90) 

70±7 

CSIMon-000441 

2008 

3.97 

0.17-0.23 (±0.11) 

160 (±20)-310 (±50) 

70±7 

CSIMon-000456 

2011 

6.92 

0.18-0.24 (±0.05) 

360 

74±3 

CSIMon-000498 

2008 

4.73 

0.17-0.21 (±0.10) 

145 (±35) - 165 (±15) 

70±6 

CSIMon-000498 

2011 

4.77 

0.19-0.22 (±0.10) 

130 (±50) - 165 (±15) 

70±6 

CSIMon-000654 

2008 

3.87 

0.10-0.13 (±0.06) 

180 (±60)-310 (±50) 

73±4 

CSIMon-000660 

2008 

7.64 

0.25-0.33 (±0.04) 

260 (±20) - 360 

74±3 

CSIMon-000660 

2011 

7.50 

0.27-0.29 (±0.05) 

200 (±30) - 360 

74±3 

CSIMon-000774 

2008 

4.73 

0.15-0.21 (±0.03) 

220 (±40) - 360 

75±2 

CSIMon-000811 

2011 

8.19 

0.26-0.28 (±0.08) 

360 

72±5 

CSIMon-000824 

2008 

7.90 

0.24-0.27 (±0.08) 

260 (±40) - 360 

72±5 

CSIMon-000928 

2008 

6.36 

0.28-0.30 (±0.10) 

120 - 300 (±60) 

71±6 

CSIMon-001054 

2011 

3.52 

0.16-0.21 (±0.09) 

270 (±70) - 360 

70±5 

CSIMon-001140 

2008 

6.79 

0.16-0.24 (±0.03) 

190 (±20) - 360 

75±2 

CSIMon-001140 

2011 

6.83 

0.14-0.24 (±0.03) 

190 (±20) - 360 

75±2 

CSIMon-001167 

2011 

5.37 

0.15-0.18 (±0.05) 

330 (±30) - 360 

74±3 

CSIMon-001296 

2011 

8.52 

0.31-0.34 (±0.13) 

260 (±80) - 330 (±30) 

69±8 

CSIMon-001308 

2008 

7.86 

0.28-0.30 (±0.06) 

330 (±30) - 360 

73±4 

CSIMon-001308 

2011 

8.04 

0.20-0.27 (±0.06) 

340 (±20) - 360 

73±4 

CSIMon-014132 

2011 

3.98 

0.18-0.32 (±0.14) 

120 (±40) - 290 (±80) 

68±9 


Notes. R w is the radius where the warp is located, assumed to be at the co-rotation radius (i.e., R w = R co )\ h max is the maximum warp height found 
for each flux dip, expressed in stellar radii R *; and (p w is the azimuthal extension of the warp found for each flux dip. The inclinations shown here 
were constrained using the occultation model. 


reproduce the widths of the observed light curves with physi¬ 
cally possible values of (f> w . This restricted the inclinations much 
more than the original calculation, resulting in much smaller er¬ 
ror bars, as can be seen in Table 3. Fig. 14 shows an example of 
the best fits for a light curve using the originally calculated value 
of inclination (with 0 W < 360°) and the new value presented in 
Table 3. It is easy to see that this value reproduces the widths of 
the flux dips much more closely. 


Warp parameters and fitted light curves The occultation 
model was used interactively to simulate the CoRoT light curves 
using the values given in Table 2 for stellar mass, radius, and co¬ 
rotation radius, with the warp maximum height and azimuthal 
extension as free parameters, and inclination as a semi-free pa¬ 
rameter, constrained to values within the error bars presented in 
Table 2. Each minimum of each light curve was fitted indepen¬ 
dently to account for slight changes in the warp from one rotation 
cycle to another. Some examples of the synthetic light curves re¬ 
produced by the occultation model are shown in Fig. 15 (syn¬ 
thetic light curves of all other stars are shown in the online ma¬ 
terial, Fig. 23). Inclinations were determined first, as described 
above, by verifying for which values we were able to reproduce 
the width of each minima without having to resort to unphysical 
values of azimuthal extension. We used a maximum inclination 
of 77°, since for inclinations larger than this the outer flared disk 


would overshadow the inner disk region 4 , making it impossible 
to observe the stellar photosphere (Bertout 2000). Error bars in 
warp maximum height and azimuthal extension were calculated 
by modeling the minima with the value of inclination fixed at 
its lowest and highest possible values, and finding the best fit in 
each case. 

Table 3 shows the final values obtained for the model pa¬ 
rameters. Since they vary from one rotation cycle to another, we 
present the intervals of the values found for each light curve. 
Most of the values found for the warp’s maximum height and 
azimuthal extension are similar to those found for AA Tau. The 
average ratio between the warp’s maximum height and the ra¬ 
dius at which it is located is h max /R w = 0.23, though individ¬ 
ual values range from 0.10 to 0.34 (for comparison we note 
that the value used for AA Tau in Bouvier et al. (1999) was of 
hmax/Rw = 0.30). The average value of the warp’s azimuthal ex¬ 
tension is 0 C = 300°, but is often a full 360°, the same as AA 
Tau, and can be as low as 110°. Inclinations generally need to 
be greater than ~ 65° for the occultation model to successfully 
reproduce the widths of the minima, but in some cases may be 
as low as 60° for relatively large warp heights and shallow flux 
dips. 


4 This effect may not be as significant for sources that have anemic 
disks, since they may possess flatter outer disks. However, this would 
only be relevant for the star Mon-314, which has aperiodic photometric 
behavior and therefore is not analyzed using the occultation model at 
any rate. 
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Fig. 15: Examples of simulated light curves of A A Tau-like stars (red) plotted over CoRoT light curves (black). 


We see that the warp maximum height should vary between 
rotation cycles by on average 11%, and at times by up to 57%, to 
account for the variations in amplitude of the light curves. The 
azimuthal extension of the warp should vary by on average 17%, 
and up to 65%, between rotation cycles to account for the varia¬ 
tions in width of the minima present in the CoRoT light curves. 
If the photometric variability of these stars is truly due to occul- 
tations by an inner disk warp caused by the interaction between 
the magnetosphere and the inner disk region, as we propose, then 
these values show how dynamic this interaction is, as has been 
predicted by MHD models. 

3.6. Color-magnitude diagrams and color variation 

We constructed color magnitude diagrams of the stars in our 
sample (both those that present AA Tau-like light curves and 
those that present aperiodic extinction dominated light curves), 
to search for trends characteristic of extinction. We used data 
taken simultaneously with CoRoT and Spitzer , plotting CoRoT 
magnitude against CoRoT magnitude minus Spitzer IRACz^m 
magnitude. It is not possible to use this method to clearly dis¬ 
tinguish if the photometric variability is due to extinction or a 
configuration of spots on the surface, but if we assume that ex¬ 
tinction is the main source of the flux dips present in both the 
CoRoT and IRAC light curves, then we can estimate extinction 
laws from the slopes of the distributions (Fig. 16). 

In order to be sure that we are analyzing parts of the light 
curve where the same phenomenon is responsible for the vari¬ 
ability in the optical and infrared bands, we calculated slopes 
using linear regression only for the parts of the data where the 
CoRoT and IRAC light curves correlate reasonably well. We also 
separated these parts in order to calculate slopes for different flux 
dips and different parts of the flux dips, as shown in Fig. 16, in 
order to verify if different extinction laws could be found for 
different parts of the disk warp. We find that entrances into and 
exits from the occultations show the same slopes, within the un¬ 
certainties. In one case, Mon-456 (see left panels of Fig. 16), 


slightly different slopes were found during different flux dips. 
This might indicate a small variability in the grain distribution 
of the inner disk warp during rotation cycles, or it could be due 
to other factors, such as variable disk emission, interfering with 
the IRAC light curves. 

We were able to estimate extinction laws for ten of the stars 
in our sample using this method, finding values between those 
consistent with ISM (A^^/Ar ~ 0.05) and up to A^ m /A R ~ 
0.29 ± 0.04. Slopes of values less than 1 were not considered, 
since extinction laws derived would be unphysical. In these 
cases, other factors, such as variable disk emission in 4.5/tm or a 
configuration of hot or cold spots on the stellar surface, could be 
strongly influencing the distribution of points in the color mag¬ 
nitude diagrams. It is possible that variable disk emission is also 
influencing our derived extinction laws, causing us to find lower 
slopes and therefore lower extinction laws than were estimated 
in Sect. 3.4. More accurate estimates will require a more exten¬ 
sive study of these and other color magnitude diagrams, which 
will be treated in a future paper. 

To check for color variation in the periodic and quasi- 
periodic stars in our sample, we compared CoRoT and I-band 
light curves with u — r for the stars whose light curves were clas¬ 
sified as AA Tau-like or AA Tau candidates in 2011. We plotted 
these light curves in phase, using the periods determined from 
the CoRoT data (Table 1). The same initial date was taken for 
all three plots for easy comparison. Though the CFHT data were 
not taken simultaneously with CoRot, by folding the light curves 
in phase, using the same initial date and period, we can see how 
the color varies during the flux dips. For most stars, we see a 
reddening at phase 0.5 during the eclipses, as is to be expected 
from an extinction event. 

For four stars (Mon-250, Mon-358, Mon-456, and Mon- 
1054), we can observe a bluing at phase 0.5 for one or two rota¬ 
tion cycles, while the other cycles either maintain relatively con¬ 
stant color or become slightly redder. This effect was observed in 
A A Tau (Bouvier et al. 2003), which showed little color variation 
except during some minima, when it became bluer. This was at- 


Article number, page 17 of 27 













































































































A&A proofs: manuscript no. AATau-like_in_NGC2264 


Mon-456 



JD - 2450000 


Mon-660 Mon-1054 



JD - 2450000 JD - 2450000 


Mon-456 


Mon-660 



14.0 


o) 14.2 
cc 


*o 144 
DC 

O 14.6 


14.8 



Slope: 1.16+- 
Slope: 1.09+- 0.02 
Slope: 1.15+- 0.03 
Slope: 1.04+- 0.02 
Slope: 1.14+- 0.02 
Slope: 1.09+- 0.03 
Slope: 1.01+- 0.05 
Slope: 1.06+- 0.02 
Slope: 1.15+- 0.05 


3.8 4.0 4.2 4.4 

CoRoT Mag - [4.5 jim] 


4.6 


Mon-1054 



Fig. 16: Top: CoRoT light curve of the stars Mon-456, Mon-660, and Mon-1054. Bottom: CoRoT vs. CoRoT-IRACz^m color 
magnitude diagrams of the same stars. The slopes were calculated for various parts of the light curves where CoRoT and Spitzer 
IRAC data best correlate. Different colors show different parts of the light curves for which slopes were calculated. 


tributed to the appearance of part of the hot spot associated with 
the accretion shocks on the stellar surface, during eclipses from 
the inner disk warp. This shows very strong evidence that these 
four stars are undergoing the same physical processes as AA Tau 
did during its stable accretion phase. 

Three stars undergo an interesting phenomenon where they 
become slightly bluer during the beginning of the eclipses and 
then become redder. This could also be evidence that the hot spot 
associated with the accretion shock is partially showing, but in 
this case it is soon occulted by the inner disk warp. These stars 
are Mon-1140, Mon-1308, and Mon-1167. For two stars, Mon- 
56 and Mon-14132, we see little or no change in color. Fig. 17 
shows some of these light curves and u — r color diagrams folded 
in phase. The others are shown in the online material, Fig. 24. 

3.7. Other characteristics of AA Tau-tike stars 

Fig. 18, top left panel, shows a mass histogram of AA Tau-like 
stars (red) and stars with aperiodic extinction dominated light 
curves (purple), plotted over a histogram of all CTTS in our 
study of NGC 2264 that present CoRoT light curves from either 
observing run. While the CTTS show a tendency toward lower 
masses, the number of stars declining rapidly for larger masses, 
the AA Tau-like stars are equally present in the whole range 
of masses considered. Therefore, though there are AA Tau-like 
stars of all masses in the CTTS range, the ratio of these stars 
to the total number of stars in that mass bin increases with in¬ 
creasing mass, meaning they are relatively more common among 
higher mass stars than lower mass ones (top right diagram in Fig. 
18). Stars that present AA Tau-like photometric behavior in ei¬ 
ther 2008 or 2011 were included in the AA Tau category. 

Inclination is a major factor in whether we observe A A Tau- 
like events or not; therefore, we know that the observed A A Tau- 
like stars are only a fraction of the total number of stars that 


present an inner disk warp. For stars of M < 0.7M o , fewer than 
10% are observed as AA Tau-like. If we consider that it is possi¬ 
ble to view the occultations at angles between ~ 60° and 77°, as 
was shown using the occultation model, then only at most 50% 
of stars in this mass range possess an inner disk warp with two 
major accretion funnels. For stars of M > 0.7M o , about 35% 
show AA Tau-like behavior, meaning that if we assume random 
inclinations, all of these stars likely possess an inner disk warp 
during some part of their evolution. 

The fraction of stars that present AA Tau-like photomet¬ 
ric behavior is quite large as we approach even higher masses 
(M > 1.0M©), much larger than we would expect if ah CTTS ob¬ 
served by CoRoT possessed an inner disk warp, given that when 
we consider random inclinations we should only see ~ 20% of 
A A Tau-like systems. This is possibly due to a selection bias, 
since the CoRoT targets are not selected randomly, but based on 
specific criteria, which included many systems with known pe¬ 
riodicity, and in 2011 included systems previously classified as 
AA Tau-like. Therefore, to compare with a less biased sample, 
we plotted a histogram of the number of AA Tau-like systems 
observed by CoRoT compared to the total number of CTTS ob¬ 
served by the CFHT MegaCam (presented in Venuti et al. 2014), 
including those that were and were not observed by CoRoT. This 
plot, shown in the bottom panel of Fig. 18, may also have a bi¬ 
ased tendency, since it is possible that some of the systems not 
observed by CoRoT also present AA Tau-like variability that 
would be unaccounted for. Even so, we see that slightly over 
20% of stars with mass M > 1.0M© present this type of variabil¬ 
ity. Extrapolating this amount to the systems where the inclina¬ 
tions prevent us from observing this tendency, we still conclude 
that nearly ah of the systems with mass 1.0M© < M < 2.0M© 
and most (~ 75%) of the systems with mass 0.7M© < M < 
1.0 M q should present an inner disk warp, while for systems of 
M < 0.7 M q this fraction would be of only ~ 25%. 
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Fig. 17: Examples of I-band and u-r light curves folded in phase. The colors used for u-r represent the same rotation cycle as the 
filled diamonds of the same color in I. Mon-56 shows no measurable variability in u-r; Mon-296 shows reddening events in u-r 
during its eclipses in the I-band; Mon-456 shows a slight reddening during one of its I-band minima (blue diamonds) and bluing 
events during the other I-band minima (orange and green diamonds); and Mon-1308 shows a bluing, followed by a reddening event 
during its minima in the CoRoT light curve. It was not possible to identify flux dips in the phase-folded I-band light curve of Mon- 
1308 because of insufficient cadence, so the phase-folded CoRoT light curve is shown instead. The same period and initial date 
were used for CoRoT and u-r, though in this case colors do not represent the same phase in both plots. 


This may be due to different stellar magnetic field configu¬ 
rations. The dipolar component of the magnetic field is respon¬ 
sible for the interaction with the inner disk, and a stronger dipo¬ 
lar component leads to a larger truncation radius. As was dis¬ 
cussed in Sect. 3.3, a larger truncation radius tends to favor a sta¬ 
ble accretion regime. Gregory et al. (2012) found that the mag¬ 
netic field configuration of a star depends strongly on its internal 
structure, and that for intermediate to high mass T Tauri Stars 
(0.5M© < M* < 2.0M©), those that are fully convective, or have 
as yet only developed a small radiative core, tend to have strong 
dipoles. At the estimated age of 2 - 3 Myr of NGC 2264, Gre¬ 


gory et al. (2012) show that stars in this mass range have not 
yet developed a significant radiative core. They should therefore 
possess strong dipolar components in their magnetic fields. 

There is still very little information available on magnetic 
fields of T Tauri stars of smaller mass (M* < 0.5M©), because of 
observational constraints. At least one T Tauri star in this mass 
range, V2247 Oph, has had its magnetic field topology measured 
and was found to have a weak dipolar component (Donati et al. 
2010). Therefore, it is possible that among fully convective stars, 
a significant difference in magnetic field configurations may ex¬ 
ist among these two mass ranges. Morin et al. (2010) showed 
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Fig. 18: Top left: mass histogram of all CTTS with CoRoT light 
curves (black), stars with A A Tau-like light curves (red), and 
stars whose CoRoT light curves are dominated by aperiodic ex¬ 
tinction events (purple). Top right: number of AA Tau-like stars 
(red) and aperiodic extinction stars (purple), divided by the to¬ 
tal number of CTTS observed by CoRoT in that mass bin. Stars 
of both CoRoT observing runs were used, and any star that pre¬ 
sented AA Tau-like photometric behavior in either run was in¬ 
cluded in the A A Tau category. Bottom: number of A A Tau-like 
stars (red) and aperiodic extinction stars (purple), divided by the 
total number of CTTS in the sample of Venuti et al. (2014), in¬ 
cluding those that were not observed by CoRoT. 


that among main sequence stars, those of very low mass show 
magnetic field topologies that differ considerably from fully con¬ 
vective stars of somewhat higher mass. While the latter generally 
possess magnetic fields with strong dipoles, main sequence stars 
of very low mass (M* < 0.2M©) show a number of different mag¬ 
netic topologies, ranging from very strong, axisymmetric, nearly 
dipolar fields to weaker fields with a strong non-axisymmetric 
component, and a strong toroidal component. This supports the 
idea that a different magnetic field regime may dominate among 
stars of lower mass. 

However, our results alone are insufficient to conclude that 
a significant difference exists among magnetic field configura¬ 
tions in these two mass ranges. There are other factors that may 
influence our higher detection of AA Tau-like systems among 
those of somewhat higher mass. For instance, mass accretion 
rates tend to decrease toward lower masses (Venuti et al. 2014), 
which would increase the truncation radius at this mass range. 
This would favor a stable accretion mechanism, but could also 
lead to the disk warp being located farther away from the co¬ 
rotation radius than we have assumed. If they are located signif¬ 
icantly farther away, we would need higher inclinations to ob¬ 
serve the eclipses in these systems, making it more difficult to 
do so. 

It is also possible that this effect is simply the result of an 
observational bias. Stars of lower mass are generally fainter than 
those of higher mass, and the signal from the CoRoT observa¬ 


tions could be too low among stars below a certain mass limit to 
allow an accurate classification of AA Tau-like light curves. 


4. Discussion 

The occultation model initially proposed for the star AA Tau ap¬ 
pears to be adequate at reproducing the general aspects of the 
flux dips of the stars listed in Table 3. It is capable of repro¬ 
ducing their widths and amplitudes using values of inclination, 
warp maximum height, and azimuthal extension that are consis¬ 
tent with those found for AA Tau itself, with only a few stars 
showing values of warp height and azimuthal extension that are 
significantly lower. 

We note that for stars whose light curves are classified as 
AA Tau-like in both CoRoT observing runs (2008 and 2011), the 
values found for warp maximum height and azimuthal extension 
in one epoch are always within the same range as those found 
in the other. This means that, although the occulting structure 
should be variable on a timescale of days to weeks because of 
the dynamic interactions, it should be stable enough to persist 
on a timescale of a few years. We also note that, for these stars, 
the periods found in each observing run are very similar to each 
other. Most are equal within their uncertainties. In all of these 
cases the radius where the occulting structure should be located 
differs by less than 4% from one epoch to the other, and most 
by less than 2%, which is within their respective uncertainties. 
It is possible that the structures move somewhat in radius over 
timescales of a few years, but they should do so by very little. 

Though it can reproduce the widths and amplitudes of light 
curve minima, the occultation model is unable to account for 
the exact shape of the flux dips, which are much more complex 
than the nearly Gaussian shapes that it reproduces. This model is 
based on a symmetric warp of homogeneous optical thickness, 
which likely does not reflect the reality of these systems. For A A 
Tau, short term brightening episodes were seen to occur during 
the occultation events, and were attributed to inhomogeneities 
of reduced optical thickness in the occulting material (Bouvier 
et al. 1999). We observe these brightening episodes in many of 
the light curves in our sample. In most cases the light curve min¬ 
ima present significant structure, which could be due to these 
inhomogeneities, to the existence of a configuration of hot spots 
visible during the occultations, or to non-steady accretion within 
the funnel flows, which would lead to variable emission during 
the flux dips. 

In some cases we observe two minima where we would ex¬ 
pect to see one, such as in the 2008 CoRoT light curve of the star 
Mon-250 (left panel of Fig. 19). We know from its 2011 CoRoT 
light curve and spectra that this star’s rotation period is on the 
order of 8.6 days (Fig. 7); therefore, the 2008 light curve spans 
only three rotation cycles, though we see more than three flux 
dips. We can speculate that the two middle dips are part of one 
structure, as the last two dips would be, and that there is signif¬ 
icant emission from the hot spot associated with the accretion 
shock superimposed onto the extinction from the base of the ac¬ 
cretion funnel, or simply that the warp was very inhomogeneous. 
Unfortunately, there is no color information simultaneous with 
this light curve to verify these scenarios. However, the occulta¬ 
tion model is incapable of reproducing one flux dip that would 
encompass both middle dips seen in the light curve using phys¬ 
ically plausible parameters. We are able to reproduce each dip 
separately, however, if we consider that there are two accretion 
funnels, separated by 160°, i.e., nearly opposite each other (right 
panel of Fig. 19). We propose that it is possible for the main ac- 
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Fig. 19: Left: 2008 CoRoT light curve of Mon-250. Right: syn¬ 
thetic light curves generated with the occultation model for two 
separate accretion funnels (blue and red) overplotted on the 2008 
CoRoT light curve of Mon-250. 

cretion funnel to at times give way to two accretion funnels as a 
result of an increase in instabilities. 

This same behavior was observed in AA Tau in 1999, when 
its quasi-periodic flux dips gave way to two flux dips per stellar 
rotation (see Fig. 15 in Bouvier et al. 2007). We may also be 
seeing this in the star Mon-1054, whose veiling features support 
the AA Tau-like scenario as the primary cause for its photomet¬ 
ric variability (see Sect. 3.2). The left panel of Fig. 20 shows its 
2011 CoRoT light curve. Arrows point to moments where we 
observe two flux dips when we would expect to see only one in a 
typical A A Tau-like scenario. The occultation model is capable 
of reproducing some of these double dips with only one accre¬ 
tion funnel, assuming that another mechanism is responsible for 
the emission within the dip, but it is not able to do so for the fifth 
double dip indicated in the left panel of Fig. 20. In this case, only 
two funnels (separated by 160° of similar height and azimuthal 
extension) can successfully reproduce both minima in the light 
curve. This may be an indication that this star is undergoing a 
transition between a stable and an unstable accretion regime, or 
that it is in an intermediate state of some instability. 

This may be the case as well for the star Mon-379 during the 
2011 CoRoT observations (Fig. 20, right panel), where we can 
observe periodic AA Tau-like dips that are easily reproduced by 
the occultation model accompanied by other extinction events. 
These additional extinction events may be due to random ac¬ 
cretion streams characteristic of an unstable accretion regime, 
as was discussed in Sect. 3.3, that coexist with the main stable 
accretion funnel. We also note that for this particular star some 
of the dips are quite narrow and the best fit found with the oc¬ 
cultation model has azimuthal extension as low as 110°, unlike 
the majority of other stars where this value is between 260° and 
360°. This may also be an indication of this star’s declining sta¬ 
bility. 

The phenomenon observed in the aperiodic light curves is 
similar to those observed in UX Ori-type objects, which exhibit 
deep, non-periodic algol-like minima caused by extinction from 
circumstellar material in the observer’s line of sight (see the 
review by Grinin 2000). Both our aperiodic stars and the UX 
Ori-type objects show evidence of being seen at high inclina¬ 
tions, which is favorable for circumstellar material located near 
the disk plane to intersect our line of sight. UX Ori flux dips 
typically present larger amplitudes and longer duration than the 
stars in our sample, probably because of larger scale disk insta¬ 
bilities. They are also much more common among Herbig AeBe 
stars, though a few CTTS have been classified as UXors (Herbst 
& Shevchenko 1999; Grinin 2000). It is possible that some of 


Fig. 20: Left: 2011 CoRoT light curve of Mon-1054. Arrows 
show moments when two flux dips are present where we would 
expect to see only one in a typical AA Tau-like scenario. Right: 
Simulated light curve (red) overplotted on the 2011 CoRoT light 
curve of Mon-379 (black). 

the stars in our sample with aperiodic behavior could be consid¬ 
ered UXors, though a study of color variability and polarization 
would be required to determine this. 

5. Conclusions 

Classical T Tauri stars that present extinction dominated light 
curves, such as the quasi-periodic A A Tau-like behavior or ape¬ 
riodic extinction events, are common in NGC 2264. Of 159 clus¬ 
ter members that show signs of accretion, 33 present one of these 
two behaviors, i.e., (21 ± 4)%. These objects give us insight into 
the phenomena that take place in the inner region of circumstel¬ 
lar disks, a region that is otherwise difficult to study because of 
its proximity to the star and to instrumental limitations. 

Extinction laws that include dust grains considerably larger 
than those found typically for the interstellar medium are neces¬ 
sary to account for the ratios between amplitudes in optical and 
IR light curves. This suggests that grain growth has occurred in 
these disks. 

A model proposed to explain the light curve of the CTTS 
AA Tau, in which at the base of an accretion column there is 
an optically thick warp in the inner circumstellar disk region 
that occults the stellar photosphere periodically, was tested on 
21 CTTS in NGC 2264. This occultation model was shown to 
be successful at reproducing the widths and amplitudes of A A 
Tau-like flux dips in the light curves of most of these stars, using 
warp parameters that are similar to those used for AA Tau. Typ¬ 
ical ratios between this warp’s maximum height and the radius 
at which it is located were found to be between h max /R w = 0.2 
and h max /R w = 0.3. The warp’s height and azimuthal extension 
were shown to be variable on a timescale of days, presenting 
typical variations of 10 - 20% between rotation cycles, though 
there is strong evidence that the occulting structure can exist over 
timescales of a few years. We have shown that, when this is the 
case, the warp maintains similar parameters during that time and 
remains at about the same radial distance. 

For the few cases where the occultation model is unable to 
reproduce an AA Tau-like light curve using a classical inner disk 
warp, we have proposed a scenario where more than one accre¬ 
tion column exists in each hemisphere, and dust is lifted above 
the disk plane at other locations around the disk, all of which oc¬ 
cult the star. This can account for the multiple flux dips per stellar 
rotation seen in the light curves. In many cases, even when a sin¬ 
gle accretion column in the visible hemisphere is able to account 
for the main features of an AA Tau-like light curve, we observe 
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various narrow, shallow dips alongside the deeper, broader flux 
dips attributed to occultation by the warp associated with the 
main accretion column. We attribute these traits to the existence 
of secondary accretion streams that coexist with the main accre¬ 
tion funnel flows. 

We ascribe AA Tau-like light curves to a stable accretion sce¬ 
nario, and light curves that present aperiodic extinction events to 
an unstable accretion scenario, and have shown that a star may 
suffer a transition from one to the other in the few years that sepa¬ 
rate our two CoRoT observing runs. For two stars with aperiodic 
extinction dominated light curves, we find evidence of hot spots 
associated with occultations, through an increase in veiling dur¬ 
ing light curve minima. We have also shown that AA Tau-like 
light curves are more common among stars of intermediate to 
high mass within the CTTS range (0.7M© < M* < 2.0Af 0 ) than 
among those of lower mass (M* < 0.7M©), an aspect we believe 
may be linked to stellar magnetic field configurations. 

We have shown that all but one star in our initial sample of 
possible AA Tau-like systems (Mon-1131, which is seen at a low 
inclination) fit very well within the A A Tau-like scenario. For 
nine of these stars we see clear evidence of a hot spot appearing 
simultaneously, or nearly so, with the occultation, either through 
an increase in veiling in the spectra or through a bluing effect 
observed in the color photometry. This strongly supports the idea 
that the flux dips observed in the CoRoT light curves are due to 
occultation by an inner disk warp that is associated with stable 
accretion funnel flows. We can therefore state with reasonable 
confidence that these are A A Tau-like systems. 
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Fig. 21: Plots of LilEW vs. I magnitude or CoRoT counts. 
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Fig. 22: Spitzer IRAC 3 .6pm and 4.5pm light curves (blue and red filled circles, respectively) overplotted on CoRoT light curves 
(black circles). The CoRoT and IRAC 3 .6pm light curves were shifted in magnitude for easier comparison. 
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Fig. 22: Continued, 
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Fig. 23: Simulated light curves of AA Tau-like stars (red) plotted over CoRoT light curves (black). 
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Fig. 24:1-band or CoRoT light curves and u-r color diagrams folded in phase. 
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Fig. 24: Continued. 

















